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Preface 



This publication describes the Fortran iv mathemati- 
cal, service, and input/output (i/o) subprograms for 
both Fortran and assembler-language programmers. 
Included is detailed information on; 

• Algorithms within the mathematical subprograms 

• Sizes of the subprograms 

• Use of the subprograms by Fortran programmers 

• Use of the mathematical and service subprograms 
by assembler language programmers 

• Techniques for replacing the tss versions of sub- 
programs with user-written versions. 
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There are also references to: 
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Introduction 



The Fortran iv library contains three types of subpro- 
grams: mathematical, service, and input/output (i/o). 
Although these subprograms are written specifically 
for Fortran programmers, they are also available to 
assembler-language programmers who use the correct 
linkage and pass the necessary information (see Ap- 
pendix B). All library subprograms are written in as- 
sembler language. 

The mathematical subprograms are similar to funo 
tion subprograms, because they are mathematical or 
computational in nature, and always return one answer 
i function value) to the calling program. Mathematical 
subprograms can be categorized by use: 

1. Direct reference, as in reference to the sine subpro- 
gram in the statement 

X = SIN(Y) 

2. Indirect reference, as in reference to an exponentia- 
tion subprogram in the statement 

X = Y**I 

The service subprograms correspond to a subpro- 
gram defined with a subroutine statement in a For- 
tran source program. These subprograms are called 
with a call statement or are implicitly called by the 
occurrence of certain situations during execution. Serv- 



ice subprograms test program-simulated machine in- 
dicators or perform utility functions. 

The Fortran i/o library consists of twenty routines 
that link together in various ways, depending upon the 
function to be performed, i/o routines are not usually 
thought of as subprograms, because any single i/o 
function depends on a number of routines. Neverthe- 
less, the Fortran i/o library can be thought of as three 
major subprograms — the Control Initialization, List 
Item Processor, and List Termination routines — -and 
seventeen supporting subprograms. This categoriza- 
tion is based upon the fact that when control is passed 
from a Fortran program to the Fortran i/o library, it 
is always one of these three routines that receives 
control. 

The Fortran i/o routines may also be categorized, 
by function, into language control routines and data 
conversion routines. These groups interact, in fulfilling 
an i/o request, by means of a communication and work 
region called the i/o Communication Routine. 

Any reference by the user program to a Fortran iv 
library subprogram causes a search of syslib for that 
program at execution time. Normally this search ob- 
tains for the user the subprogram provided as part of 
tss. The user can, however, provide his own version 
of the subprogram, as described in Appendix A. 
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Section 1: Mathematical Subprograms 



The two types of mathematical subprograms are di- 
rectly referenced subprograms and indirectly refer- 
enced subprograms. The directly referenced subpro- 
grams are called by the object program in response to 
a statement of the form 

X=SIN (Y) 

In this statement, direct reference is made to the 
mathematical sine subprogram, by its entry name; sin. 

An example of indirectly referenced subprogram 
usage 1 is the call on an exponentiation subprogram, 
made as the result of a statement of the form 
X = Y * * I 

In this statement, no direct reference is made to a 
subprogram by the Fortran programmer. The Fortran 
compiler determines that a subprogram is required to 
perform the exponentiation operation, however, and 
causes the object program to call the appropriate ex- 
ponentiation subprogram. 

The algorithms describing the method of computa- 
tion of the* mathematical subprograms are given in 
Appendix F. Other information concerning these sub- 
programs is contained in Tables 1, 2, 3, 4, and 5 of this 
section and in Appendix A. 

Tables 1 and 2 give this information: 

Function: A brief description of the type of mathe- 
matical operation periormed. 

Entrif Name: The mathematical subprograms contain 
an entry point corresponding to each name that may be 
directly referenced (such as sin) and each name that 
may be indirectly referenced (such as ciicbga, when 
raising an 1*4 integer to an I * 4 power). This column 
shows all entry points in the mathematical subpro- 
grams. 

Definition: This column gives a mathematical equa- 
tion that represents the computation. (It is not meant 
to represent the way the subprogram is called.) An 
alternative equation is given when there is another way 
of representing the computation in mathematical nota- 
tion. For example, the square root can be represented 
as either 

Vx or x 1/2 

Arguments): These columns describe the value(s) 
for which the function value is to be computed. 

* Argument Number — The number of arguments (one 
or two ) that the user must supply. 

• Argument Type — The type and length of each argu- 
ment. Integer, real, and complex represent the type 



of number; the notations *4, *8, and *16 represent 
the length, in bytes, of the argument. 
Note: In Fortran iv, a real argument corresponds 
to the real°4 argument, and a double-precision argu- 
ment corresponds to the real* 8 argument. A single- 
precision complex argument corresponds to the com- 
plex* 8 argument, and a double-precision complex argu- 
ment corresponds to the complex *i6^rgument. 
• Argument Range (Table 1 only) — The valid range 
for each argument. If an argument is not within its 
valid range, an error message is issued and execution 
of this load module is terminated. (See the Error 
Condition and Error Message column descriptions 
below. ) 

Function Value Returned: This column describes the 
function value returned by the subprogram; the nota- 
tion is the same as that used for the argument type. 

Error Condition: This column describes the argu- 
ment ranges not allowed when using the mathematical 
subprogram. 

Storage Estimates: This column shows the approxi- 
mate number of bytes required for each mathematical 
subprogram; the approximate, total size of each sub- 
program's csect and psect. ( Fortran iv mathematical 
subprograms each contain one public, read-only, re- 
enterable csect and one psect. The length of each of 
the control sections is less than 4096 bytes. The subpro- 
grams are link edited, and their csects are combined.) 

Other Subprograms Required: Many mathematical 
subprograms require other mathematical subprograms 
to perform their function. The entry names of the other 
subprograms are listed in this column. (This column 
does not include chcbza, which is called by all mathe- 
matical subprograms where error exit is possible. ) 

Routine Name: Each mathematical subprogram is 
assigned a routine name that is normally of no interest 
to Fortran programmers. Appendixes A and B describe 
use of this name. 

Accuracy Figures (Table 1 only): These columns 
give accuracy figures for one or more representative 
segments within the valid argument range. The accu- 
racy figures are based upon the assumption that the 
arguments are perfect (i.e., without error and, there- 
fore, having no error-propagation effect on the an- 
swers). The only errors in the answers are those 
introduced by the subprograms. Information given in 
the accuracy-figures columns is: 



Table 1. Summary of Directly Referenced Mathematical Subprograms 
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Entry 
Name 


Definition 


Argument(j) 


Function 

Value 

Returned 


Error Condition 


Storage 
Estimates 


Other Subprograms Required 


Routine 
Nome 


Accuracy Figure! 


Argument 
Range 


Sample 

EAJ 


relat>"* 






No, 


Type 


Range 


Hex 


Dec 


M(«) 


• («) 


M(f) 


»<8) 


COMMON AND 

NATURAL 

LOGARITHM 


CDLOG 


Ln (arg) or Log e (arg) 
See Note 8 


1 


COMPLEX * 16 


org * + Oi 


COMPLEX • 16 


Argument "0 + 01 


1E8 


488 


CDABS, DLOG, DATAN2, DSQRT 


CHCAP 


The full range 
except (1 + Oi) 


Note 1 


2.72 x \0" ]6 


5.38x IO* 17 






CLOG 


Ln (org) or Log e (arg) 

See Note 8 


1 


COMPLEX * 8 


arg * + 01 


COMPLEX * 8 


Argument -0 + 01 


1D0 


464 


CABS, ALOG, ATAN2, SORT 


CHCAO 


The full ronge 
except (1 + Oi) 


Note 1 


7.15x 10" 7 


l,36x 10" 7 






DLOG 


Ln (arg) or Log e (org) 


1 
1 


REAL * 8 


arg > 


REAL * 8 


Argument 5 


21A 


538 




CHCAF 


0,5 5 X S 1.5 


U 






4.60 x IO" 17 


2.09x IO"' 7 


X-cO.S, X > 1.5 


E 


3.32x 10" 16 


5.52x IO*' 7 






DLOG10 


Log 10 (arg) 


REAL * 8 


arg > 


REAL * 8 


Argument £ 


21A 


538 




CHCAF 


0.5 S X* 1.5 


U 






2.73 k JO" 17 


1.07k IO"' 7 


X<0.5, X > 1.5 


E 


3.02x 10" 16 


6.65 x IO* 17 








ALOG 


Ln (org) or Log e (org) 


1 


REAL * 4 


arg > 


REAL * 4 


Argument S 


1D0 


464 




CHCAE 


0.5 S X $1.5 


U 






6.85 x 10" 8 


2.33 x IO* 8 


X< 0.5, X >1.5 


E 


8-32 x 10* 7 


1.19x 10" 7 












J 


arg > 






1D0 


464 




CHCAE 


0.5 i X Si. 5 


U 






7 13 x IO* 8 


2.26 x IO* 8 


ALOGIO | Log, (org) 1 i REAL • 4 

1 A 1 i. 


REAL * 4 Argument < 


X< 0.5, X > 1.5 


E 


1.05 x IO* 6 


2.17x 10*' 






ICDEXP | e arg | 1 j COMPLEX * 16 


real arg < 174.673 
jimag arg| < I w 


COMPLEX * 16 Real Argument > 174,673 

(imaginary Argument| > 2^ w 


270 


624 


DEXP, DSIN, DCOS 


CHCAN 


|Xl|$MX2|Sf 


U 


3.76x 10" 16 


l.lOx io" 16 








! 


|Xl| ^ 20, |X2)< 20 


u 


2,74x 10* 15 


9.64x 10' 16 








CEXP | e arg 1 | COMPLEX » 8 real arg 174.673 


COMPLEX * 8 iReal Argument > 174,673 

((imaginary Argumentj > 2 18r 


250 


592 


EXP, SIN, COS 


CHCAM 


N<,7o,N<.f 


u 


9.93 x I0" 7 


2.67 x 10" 7 






EXPONENTIAL 


I I 
1 I 1 


I. i ^ 8 
prnag argj < I tr 


|X,|S170, 

! < ! x §l*2o 


u 


1.07* 10~ 6 


2.73 x 10" 7 






DEXP e arg 1 1 REAL * 8 

■ 


arg < 174.673 


REAL * 8 '[Argument > 174.673 2C0 


704 




CHCAD 


|X|Sl 


u 


2.04x 10" 16 


5.43x 10" 17 












1 < | X| 5 20 


u 


2. 03« 10* ,A 


4.87x 10"' 7 














20 < | X | S 1 70 


u 


1.97x l(f ' 6 


4.98x 10'' 7 






EXP 


e ar 9 


1 


REAL * 4 1 arg < 174.673 

J 


REAL * 4 


Argument > 174.673 


1A8 424 





CHCAC 


|x|<i 


u 


4.65 x 10" 7 


1.28 x 10" 7 






1< |X| |170 


u 


4,42 x 10~ 7 


1.15 x 10" 7 






SQUARE ROOT 


CDSQRT 


1 2 / 
(org) or ^J org 


1 


COMPLEX * 16 


Any 


COMPLEX * 16 


None 


148 


328 CDABS, DSQRT 


CHCAT 


The full range 


Note 1 


1.76x IO* 16 


4.06x )0" 17 






CSQRT 


(org) or +J arg 


1 


COMPLEX * 8 


Any 


COMPLEX * 8 


None 


138 


312 CABS, SORT 


CHCAS 


The full range 


Note 1 


7.00x 10" 7 


1.71 x 10" 7 






DSQRT 


(arg) ' or *J arg 


1 


REAL * 8 


arg # 


REAL * 8 


Negative Argument 


160 


352 




CHCAB 


The full range 


E 


1.06x 10" 16 


2.16x 10" 17 






SORT 


(arg) or *J org 


1 


REAL * 4 


org #0 


REAL * 4 


Negative Argument 


158 


344 




CHCAA 


The full range 


E 


4.45 x 10" 7 


8.43 x 10" 8 






ARCSINE 

AND 

ARCCOSINE 


DARSIN 


arcsine (arg) 


1 


REAL * 8 


|arg| < 1 


REAL * 8 


Argument > 1 


288 


648 


DSQRT 


CHCAX 


-1 < X <+l 


u 


2.04 x 10" 16 


5.15x 10" 17 






DARCOS 


arcco$!ne (org) 


1 


REAL * 8 


lorgl <1 


REAL * 8 


Argument > 1 


288 


648 


DSQRT 


CHCAX 


-i < x S+i 


U 


2.07x l<f' 6 


7.05 x 10" 17 






ARSIN 


arcsine (org) 


1 


REAL * 4 


arg <1 


REAL * 4 


Argument > 1 


1F0 


496 


• SORT 


CHCAW 


-i < x S + i 


U 


9.34x 10* 7 


2.06x 10" 7 






ARCOS 


arccosine (org) 


1 


REAL * 4 


|arg| <1 


REAL * 4 


|Argument| > 1 


1F0 


496 


SQRT 


CHCAW 


-1 < X < + l 


U 


8.85 x 10' 7 


3.19x 10" 7 






ARCTANGENT 


DATAN 


arcran (arg) 


1 


REAL * 8 


Any 


REAL * 8 

i _ 


None 


288 


648 




CHCBR 


The full range 


Note 7 


2.18 k 10~ 16 


7.04x 10* 17 






DATAN2 


arctan (arg| arg2) 


2 


REAL * 8 


arg # REAL * 8 


x, = x 2 -- o 


288 


648 




CHCBR 


The full range 


Note 7 


2.18* 10' 16 


7.04x 10"' 7 






ATAN 


arctan (org) 


1 


REAL * 4 


Any | REAL * 4 None 


1E8 


488 




CHCBQ 


The full range 


Note 7 


1.01 x I0' 6 


4.68x 10' 7 






ATAN2 


arctan (arg^ ''arg-) 


2 


REAL ' 4 


arg *0 


REAL * 4 


x, = x 2 -0 


1E8 


488 




CHCBQ 


The full range 


Note 7 


1.01 x io" 6 


4,68 x 10* 7 






TRIGONOMETRIC 
SINE & COSINE 


CDSIN 


iin (arg), arg in 
.adians 


1 

1 

I 


COMPLEX * 16 


imag arg < 174,673 


COMPLEX * 16 


1 |Rcal Argument] > 2^ w 

Imaginary Argument > 174,673 


340 


832 


DSIN, DCOS, DEXP 


CHCAR 


|x,|Su>,|x 2 |Si 


U 


2.35 x 10"' 5 
See Note 4 


2.25 x 10' 16 






CDCOS 


cos (arg) , arg in 
radians 


COMPLEX * 16 


|, M l.,g| <2 5 "' 

imag aig < 1 74 .673 


COMPLEX * 16 


|Real Argument| * 2^ * 
Imaginary Argument! > 174,673 


340 


832 


DSIN, DCOS, DEXP 


CHCAR 


|X,| i io,|x 2 |5i 


U 


3.98 x 10" 15 
See Note 3 


2.50 x 10"' 6 






CSIN 


sin (arg) , arg in 


1 


COMPLEX * 8 


|,ealarg| < 2 18 * 

1 imag argj < 174.673 


COMPLEX * 8 


|Real Argument > ? ,8 » 
llmaginary Atgument > 174,673 


2F8 


760 


SIN, COS, EXP 


CHCAQ 


|X,| i 10,|X^1 


U 


1.92 x IO" 6 
See Note 6 


7.38 x 10~ 7 
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Table 1. Summary of Directly Referenced Mathematical Subprograms (cont.) 

1 ? 3 4 





Fnf.w 








Argument!*) 


Function 




Storo 


3e 




I 






Accuracy Figures 


Function " ' Definition 

Nome ). 

No. 


' Value fcnoi Condition 
\ ' - -— • R efu , ned 
Type Range 


~ it '™ 1 "" m Other Subprog-amj Required 
Hex ' Dec 


Name [ Argument 
[ Range 


Sample 
E/U 


Relative 
M («) 9 ( . ) 


Absolute 
M(E) j »(E) 


CCOS 


CO* iQrg), org in 1 
rodiani 


COMPLEX * 8 


1 
|.eal org| < 2 I8 ~ ; COMPLEX * 8 

i mag org < 174,673 i 


1 Real Argument I >2 18 * 
Imaginary Argument > 174,673 


2F8 

1 


760 SIN, COS, EXP 


CHCAQ 


x,| S»o,|x 2 |< l 


U 


2.50 x IO" 6 
See Note 2 


7.66x 10** 7 






DSIN 


iin (org), org in , 1 
rodiani I 


REAL * 8 


1 1 <->50» 
jorg < 7 


REAL * 8 


Urgumentl > 2 50 " 


2B8 


696 




CHCAJ 


l*|*j 


U 


3.60k IO" 16 


4.82 x IO" 7 


7.74.10- 17 


1.98x 10* 17 


j 


i<|x|Slo 


U 






l.Mx IO* 16 


6.49 x IO" 17 








10 < | x| < 100 


U 






2.68 x IO" 15 


1.03x 10* 15 






1 


REAL * 8 


|org| < 2 50 » 


REAL * 8 


Argument 2 2 w 


288 


696 




CHCAJ 


i x< * 


U 






1.79x IO* 16 


6.53x IO** 17 


TRIGONOMETRIC ~~~~ !—'-*'< ~« - 
SINE 4 COSINE jrcxi.ons 


- 10 S x < 
w < X < 10 


U 






1.75 x IO* 16 


5.93x IO"' 7 




i 


10 c | x| £ 100 


U 






2.64 x IO* 15 


1.01 x io" 15 




SIN | lin (org), org in 
■ rodiofis 


1 


REAL * 4 


|org| <2 18 » 


REAL • 4 


|Argument| > 2 18 * 


1F8 


504 




CHCAI 


x|S| 


U 


1.32 x IO" 6 


1.82 x 10** 7 


1.18 x 10" 7 


4.55 x 10" 8 


. j 


1 <| xl i 10 

2 ' ' 


u 






1.15 x IO* 7 


4.64 x IO" 8 






10 < | x | S 100 


u 






l.Mx IO* 7 


4.52 x IO" 8 


COS 


COS (org), org in 
rodiani 


1 


REAL * 4 


|org| <2 18 ' 


REAL * 4 


Urgumentl > 2 18 * 


1F8 


504 




CHCAI 


o<.xS 


u 






1.19x IO" 7 


4.60 x IO" 8 


-io £ x < o 

w < x £ io 


u 






l.Mx IO' 7 


4.55 x IO* 8 


.. 




10 < | X | Sioo 


u 






1.14 x 10" 7 


4.60 x 10" 8 




DTAN 


tan (org), org in 
radians 


1 


REAL # 8 


|org| < 2 50 ' 


REAL * 8 


| Argument | > 2 50 * 

Argument too close to a Singularity 
(i.e. , too close to an odd 
multiple of tt/2) 


2F8 


760 




CHCAZ 


*I S J 


u 


3.41k IO* 16 


6.27 x 10" 17 






i 


?<l*l s f 


u 


-12 
1.43x 10 '* 

See Note 5 


2.95 x IO** 14 






f <|x|-io 


u 


2.78 x 10~ 13 
See Note 5 


7.23 x IO* 15 










10 <| x| JlOO 


u 


3.79 x 10" 12 
See Note 5 


9.50x 10* 14 








OCOTAN 


cotan (org), arg in 
rodians 


1 


REAL * 8 


|.rg| - 7™* 


REAL * 8 


1 1 50 
Argument j ; 2 * 

Argument too close to a Singularity 
(i.e. , too close to a 
multiple of tr) 


2F8 


760 




CHCAZ 


1*1* J 


u 


2.46x IO" 16 
See Note 5 


8.79x IO"* 17 






TRIGONOMETRIC 
TANGENT 




f<l x l s f 


u 


2.78x 10 IJ 
See Note 5 


8.61 x 10" 15 






f<|x|S 10 


u 


-l3 
5.40x 10 ,J 

See Note 5 


i 
1.13x IO* 14 


- ' ~ "' 




10 < | x| i 100 


u 


-13 
8.61 x 10 4 

See Note 5 


4.61 x IO*" 14 






TAN 


tan (org), arg in 
radiant 


1 


REAL * 4 


|org| < 2 18 * 


REAL * 4 


Urgumentl > 2 18 * 

Argument too close to a Singularity 
(i.e. , too close to an odd 
multiple of tt/2) 


288 


648 




CHCAY 


!*l* 4 ' 


u 


1.71 x IO" 6 


2.64 x IO*" 7 






* < 1 v I < » 


u 


1.05 x IO" 6 
See Note 5 


3.59 x 10" 7 






J <| x\i 10 


u 


6.49 x IO** 6 
See Note 5 


3.38 x IO* 7 






10<| X | i 100 


u 


1.57k IO** 6 
See Note 5 


3.07 x IO" 7 






COTAN 


cotan (arg), arg in 
rodians 


1 


REAL • 4 


|org| < 2 18 ' 


REAL * 4 


Argumentj > 2 ]Bw 

Argument too close to a Singularity 
(i.e. , too close to a 
multiple of tt) 




648 




CHCAY 


1*1*} 


u 


1.07 x IO* 6 


3.58 x IO* 7 






J<|x|S 5 


u 


1.40x IO* 6 
See Note 5 


2.56 x IO* 7 






£<|x|Sl0 


u 


1.30* IO" 6 

See Note 5 


3.11 x IO** 7 






10<| x| - 100 


u 


1.49 x 10" 6 
See Note 5 


3.15x 10" 7 






HYPERBOLIC 
SINE & COSINE 


DSINH 


sinh (org) 


1 


REAL * 8 


|arg|< 175,366 


REAL * 8 


| Argument 1 > 174,673 


250 


592 


DEXP 


CHCBB 


|x| i 0.88137 


u 


2.06 x IO* 16 


3.74 x IO*"' 7 






0.88137 <| X | S 5 


u 


3.80 x IO* 16 


9.21 x 10** 17 






DCOSH 


cosh (org) 


1 


REAL * 8 


|arg|< 175.366 


REAL * 8 


Argument t 174,673 


250 


592 


DEXP 


CHCBB 


-5 S X < +5 


u 


3.63 x 10* 16 


9.05 x 10" 17 






SINH 


sinh (arg) 


1 


REAL • 4 


|arg| < 175.366 


REAL * 4 


Argument i 174,673 


1F8 


504 


EXP 


CHCBA 


-5 £ X < + 5 


u 


1.26 x IO* 6 


2.17x 10" 7 






COSH 


cosh (arg) 


1 


REAL * 4 


|org|< 175.366 


REAL * 4 


Urgumentl > 174,673 


1F8 


504 


EXP 


CHCBA 


-5 < X < +5 


u 


1.27 x IO" 6 


2.63 x 10* 







Table 1. Summary of Directly Referenced Mathematical Subprograms (cont.) 



1 


2 


3 






4 




3 


6 


7 




8 


9 








10 






Function 


Entry 
Name 


Definition 


Type 


igumentiii 

Range 


function | 

Value Error Condition 

Returned 

; 


Storage 

Estimates 


r\i\. c l d • _i Routine 

Other Subprograms Required 

Name 


Accuracy Figures 


Argument 
Range 


Sample 
E/U 




absolute 


Hex 


Dec 


M (f ) 


*(0 


M(E) 


*!t) 


HYPERBOLIC 

TANGENT 


DTANH 


tonh (org) 1 


REAL * 8 


Any 

Any 

Any 

See Note 9 




130 


304 


DEXP CHCAL 


|x| i 0,54931 


U [ 1.91 x 10' 16 


3.86 x 10" 17 










0,54931 < |x| i 5 j 
MS 0.7 


U 1 .54 x 10 


-17 
1.87 k 10 






TANH 


tanh (arg) 1 i 


REAL * 4 


1A4 


•\Kh 


FXP i 


CHCAK > 


U 8.48 x 10" 7 


1.48* 10" 7 












0.7 < |x| < 5 


U ! 2.44 x !0~ 7 


4.23* 10" 8 






ABSOLUTE 
VALUE 


CDABS 


II ! 
a.g 1 


COMPLEX * 16 




REAL * 8 


C8 200 DSQRT CHCAV 


The full range 


Note 1 


2.03* 10" 16 


4.83x 10' 17 






CABS 


arg| 




COMPLEX * 8 


Any 

See Note 9 




REAL * 4 j 


CO 


192 


SORT ! CHCAU 


The full range 


Note 1 


9,15x 10~ 7 


2.00x 10" 7 






ERROR 

FUNCTION 


ERF 

DERF 


r* -J 


REAL * 4 


Any 




REAL * 4 1 


208 


520 


EXP i CHCBU 
1 


M < 1 j 


U 


8.16k 10' 7 


1.10 k io' 7 






Jfc h 


. 


1 < |x| < 2.04 


U 


1.13k 10~ 7 


3.70x 10* 8 






?.04< !xj< 3.9192 


u 


5,95 x 10" 8 


3.41 x 10 ° 






rX 2 


Real * 8 


Any 


| Real * 8 


328 808 


DEXP ! CHCBW 


|x| < 1 ! U 


189* 10" 16 1 2-60x 10' 17 






2 f e" U du 


1 


l< |x|< 2.04 


u 


2.87 x 10" 17 9.84x 10' 18 








| 


2.04 <|x| < 6.0<?2 


U j l,39x 10" 17 8.02x I0" 18 






COMPLEMENTED 
ERROR FUNCTION 


ERFC . l-erf(x)or 


] 


REAL * 4 


Any 


REAL * 4 j 

1 


208 j 520 


EXP CHCBU 


-3.8 < x < U 9 ,Q, in" 7 7 96 « '0' 7 






i X S 1 U 7.42x 10" 7 I 1.27 x I0' 7 










1 < xS 2.04 u 1 54 x 10" 6 , 3.78 * 10" 






j 2.04 < X < 4 U 
4 < X < 13.3 U 


2.28 x 10" 6 j 3.70 x 10" 






c ' f " ■; 

1.55 x 10 8.57x 10 






OERFC 


1 - erf (x) or 

I f <Vdu 1 
JrJ* € ! 


Real * 8 


Any 




Real * 8 j 


1 ' ' 


-6 < X< 


o 2 08 x ID" 16 j 6.52 x 10" 17 






j 




i X i 1 


U 1.40x lO"' 6 i 2 59x !0" 17 






328 


808 


DEXP 


CHCBW 


1 < X < 2.04 


U 1 4.11 x 10" 16 | 8.86x 10" 17 






2.04 < X<4 


U 3.26x 10" 16 : 8.65x lO* 17 










4 < X < 13.3 


U 3.51 x 10" 15 1 1.96 x 10" 15 






GAMMA (r) 


GAMMA 


f 90 u-'e-du 


REAL * 4 


X> 2" 252 
X < 57.5744 




350 i R4ft 


EXP, ALOG 


CHCBT 


< X < 1 


U 


-7 
9 86 x 10 


3.66 x 10* 7 






I u e du 
J Q 




Real Argument < 2~ 252 

j 






1 < X < 2 


U 


1.13m 10* 7 


3. 22* 10' 8 










2 < X < 4 


U 


9,47 x 10** 7 


3.79 x 10* 7 






4 < X < 8 


U 


2.26x 10" 6 


8.32* 10" 7 






8 < X < 16 


U 


2.20 k 10' 5 


7.61 * 10" 6 






16 < X < 57 


U 


4 62« 10" 5 


1.51 * 10* 5 






DGAMMA 


/ u x ' e du ; 1 

o 


REAL # 8 


X > 2~ 25? 
X < 57.5744 






420 


1056 


DEXP, DLOG 


CHCBV 


< X < 1 


U 


2 Mm 10" 16 


7.84* 10* 17 






REAL * 8 (Real Argument > 57,5744 
JReal Argument < 2' 252 

I 


1 i X < 2 


U 


h 2.52x,0" 17 


6.07* 10' 18 






2 < X < 4 


U 


2 21* 10' 16 


8.49* 10' 17 






4 < X< 8 


u 


5.05* 10"' 6 


1.90x 10 






- 






1 


8 < X < 16 


U 


6.02* 10" 15 


1.78x 10* ]5 






16 < X < 57 


U 


1.16* 10" U 


4.11 * 10" 15 






LOG - GAMMA 


ALGAMA 


r oo x -l -u 




REAL * 4 


X > 

X < 4.29)3 x 10 73 

X > 

X < 4,2913 x I0 73 




REAL * 4 (Real Argument > 4.2937 x 10 73 
Real Argument < 


350 


848 


EXP, ALOG 


CHCBT 


< X < 0.5 


U 


1.16* 10" 6 


3.54 m 10" 7 






log / u e du 1 

1 

. __|_ 


0.5 < X < 3 


U 






9.43 x 10* 7 


3.42 x !0* 7 


3 < X < 8 


U 


1.25* 10" 6 


3.04 k 10' 7 






8 i X < 16 


U 


1.18x 10* 6 


3.80* 10* 7 








| 


16 < X < 500 


u 


9.85 x 10* 


1.90* 10* 7 






DLGAMA 


log/ u e du ; 

'; 

j 


REAL * 8 


OTA, * O B»nl imnm.nt "5 A ?V%7 w If/ 


420 


1056 


DEXP, DLOG 


CHCBV 


< X i 0.5 


U 


2.77* 10' 16 


9.75* IO*' 7 








0.5 < X < 3 


U 






2.24 x 10" 16 


7.77* 10' 17 








3 i X < 8 


U 


2.89x 10" 16 


8.80* 10' 17 






8 £ X < 16 


U 


2.86x 10* 16 


8.92* 10"' 7 






16 < X < 500 


U 


-16 
1 . 99 * 10 


3.93x IO"' 7 








Notes 1 
2 


The distribution of sample org 
these statistics ore bosed is e 
and is uniform around the ori 

The maximum relative error c 
function is based upon a set c 
arguments within the range, 
proximity of the points (n + 1 
n = 0, ±1, ±2, .. .,) the relot 
quite high, although the absc 


uments upon which 
uponentiol radially 
Jin. 

ited for the CCOS 
f 2000 random 
In the immediate 
/2)tr +0i (where 
ive error can be 
lute error is ynall, 


3 
4 


The maximum relative error cited for the CDCOS 
function is based upon a set of 1500 random arguments 
within the range. In the immediate proximity of the 
points (n + 1/2)* + Oi (where n = 0, ±1, ±2, ...,) 
the relative error can be quite high although the 
absolute error is small. 

The maximum relative error cited for the CDSIN 
function is based upon a set of 1500 random arguments 
within the range. In the immediate proximity of the 
points ni + Oi (where n •= t 1, i 2, , , , ,) the relative 
error can be quite high although the absolute error is small. 


5, The figures cited as the maximum relative errors are 
those encountered in a sample of 2500 random arguments 
within the respective ranges. See the appropriate 
section in Appendix F for a description of the behnvioi 
of errors when the argument is near a singularity or o 
zero of the function. 

6. The maximum relative error cited for the C5IN function 
is based upon a set of 2000 random arguments within the 
range. In the immediate proximity of the points 

n w * Oi (where n = i 1 , ± 2, , . , ,) the relative error 
can be quite high although the absolute error is small. 


7. The somple orguments were tangents of numbers 
uniformly distributed between - w/2 end ♦ »/2. 

8. The answer given is th* principal value ,!.«., 
the one whos# imaginary port lits btnw#tn 

- tr and + tr . 

9. Floating-point overflow can occur. 
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Table 2. Summary of Indirect!}' Referenced Mathematical Subprograms 





2 


2 




4 


5 


6 


7 




8 


10 


Funcfion 


Entry 
Name 


Definition 


Argument(s) 


Function 

Value 

Returned 


Error Condition 


Storage 
Estimates 


Other Subprograms Required 


Routine 
Name 


No, 


Type 


Hex 


Dec 


RAISE AN INTEGER 
BASE TO AN [ 
INTEGER POWER 


CHCBGA 


y = i * * J 


2 

| 


i - INTEGER * 4 


INTEGER * 4 


Base !s zero 


1B4 


436 




CHCBG 




CHCBGB 


y = i * * i 


2 ! 


j ^ INTEGER * 2 


INTEGER * 2 


Base is zero and exponent is 
zero or negative 


1B4 


436 




CHCBG 




CHCBGC 


y = j * * ; 


2 I 


j = INTEGER * 2 
i = INTEGER * 4 


INTEGER * 4 


Base is zero and exponent is 
zero or negative 


1B4 


436 




CHCBG 




CHCBGD 


y - i * * j 


2 , 


i - INTEGER *T| 

j - INTEGER * 2 

1 


INTEGER * 4 


Base is zero and exponent is 
zero or negative 


1B4 


436 




CHCBG 




RAISE A REAL 

BASE TO AN 
INTEGER POWER ! 


CHCBHA 


y = a " * i 


2 ' 


a REAL * 4 j 
i - INTEGER * 4 j 


REAL * 4 


Base is zero and exponent is 
zero or negative 


144 


324 




CHCBH 




CHCBHB 


y = a * * j 


2 


a = REAL * 4 REAL * 4 

j = INTEGER * 2 

1 ., 


Base Is zero and exponent is 
zero or negative 


144 


324 




CHCBH 




RAISE A DOUBLE \ 
PRECISION BASE TO j 
AN INTEGER POWER ; 


CHCBIA 


y = o * * i 


2 ; 


a - REAL * 8 

- INTEGER * 4 


REAL * 8 


Base is zero and exponent Is 

zero or negative 


14C 


332 




CHCBI 




CHCBIB 


y " a * * j 


2 : 


a = REAL * 8 
j = INTEGER * 2 


REAL * 8 


Base is zero and exponent is 
zero or negative 


14C 


332 




CHCB1 




RAISE A REAL 
BASE TO A REAL 
POWER 


CHCBJA 


y = a * * b 


2 ' 


a = REAL * 4 
b = REAL * 4 


REAL * 4 


Base is zero and exponent is 
zero or negative 


ICO 


448 


EXP ALOG 


CHCBJ 


RAISE AN INTEGER 
BASE TO A REAL 

POWER 


CHCBJB 


y = i ** b 


2 


b = REAL * 4 
i = INTEGER * 2 


REAL * 4 


Base is zero and exponent is 
zero or negative 


ICO 


448 


EXP, ALOG 


CHCBJ 


CHCB.IC 


y = i **b 


2 


b = REAL * 4 
i = INTEGER * 4 


REAL * 4 


Base is zero and exponent is 
zero or negative 


ICO 


448 


EXP, ALOG 


CHCBJ 


RAISE A REAL OR 
INTEGER BASE TO 
A REAL POWER; 
BASE AND/OR 
EXPONENT 
DOUBLE 
PRECISION 


CHCBKA 


y =a * * b 


2 


a = REAL * 8 
b = REAL * 8 


REAL ' 8 


Base is zero and exponent is 
zero or negative 


230 


560 


DEXP, DLOG 


CHCBK 


CHCBKB 


y = j* *b 


2 


b = REAL * 8 
i = INTEGER * 2 


REAL * 8 


Base is zero and exponent is 
zero or negative 


230 


560 


DEXD, DLOG 


CHCBK 


CHCBKC 


y = i * * b 


2 


b = REAL * 8 
i = INTEGER * 4 


REAL * 8 


Base is zero and exponent is 

zero or negative 


230 


560 


DEXP, DLOG 


CHCBK 


CHCBKD 


y = a * * b 


_ 1 


a * REAL * 4 
b = REAL * 8 


REAL * 8 

See Note, 


Base is zero and exponent is 
zero or negative 


230 


560 


DEXP, DLOG 


CHCBK 


CHCBKE 


y =a * * b 


2 


a - REAL * 8 
b = REAL * 4 


REAL * 8 


Base Is zero and exponent is 
zero or negative 


230 


560 


DEXP, DLOG 


CHCBK 


RAISE A COMPLEX 
BASE TO AN 
INTEGER POWER 


CHCBMA 


y =a * * i 


2 


a = COMPLEX * 16 
i = INTEGER * 4 


COMPLEX * 16 


Base is zero and exponent is 
zero or negative 


274 


628 




CHCBM 




CHCBMB 


y = a * * j 


2 


a = COMPLEX * 16 
I = INTEGER * 2 


COMPLEX * 16 


Base is zero and exponent is 
zero or negative 


274 


628 




CHCBM 




CHCBCA 


y = a * * \ 


2 


a = COMPLEX * 8 
i = INTEGER * 4 


COMPLEX * 8 


Base is zero and exponent is 
zero or negative 


24C 


588 




CHCBC 




CHCBCB 


y = a * * j 


2 


a -COMPLEX * 8 
i = INTEGER * 2 


COMPLEX * 8 


Base is zero and exponent is 
zero or negative 


24C 


588 




CHCBC 




PRODUCE ERROR 
MESSAGE AND 
TERMINATE 
EXECUTION 


CHCBZA 












E8 


232 


As required by use of the 
EXIT macro Instruction 


CHCBZ 



NOTE: The REAL # 8 function value returned by CHCBKD is not more accurate than the REALM base 
given as an argument. 



• Argument range — This column gives the argument 
range used to obtain the accuracy figures. For each 
function, accuracy figures are given more represen- 
tative segments within the valid argument range. 
These figures are the most meaningful to the func- 
tion and range under consideration. 

The maximum relative error and standard deviation 
of the relative error are generally useful and reveal- 
ing statistics. However, they are useless for the range < 
of a function where its value becomes 0, because the 
slightest error in the argument can cause an unpre- 
dictable fluctuation in the magnitude of the answer. 
When a small argument error would have this effect, 
the maximum absolute error and standard deviation 
of the absolute error are given for the range. For ex- 
ample, absolute error is given for sin (x) for values 
of x near -w. 

• Sample — This column indicates the type of sample 
used for the accuracy figures; the type depends upon 
the function and range under consideration. The sta- 
tistics may be based either upon an exponentially 
distributed ( E ) argument sample or a uniformly dis- 
tributed ( U ) argument sample. 



Error Message: chcbzioo is issued each time an error 
occurs. This message gives the error condition, the 
entry name, and the address of the call to the math 
routine in the users program. 



Table 3, Exponentiation With Integer Base and Exponent 



Base (I) 




Exponent (J) 




J>o 


J = o 


J<o 


I>1 


Compute the 
function value 


Function 
value = 1 


Function 
value = 


1 = 1 


Compute the 
function value 


Function 
value = 1 


Function 
value = 1 


1 = 


Function 
value = 


Error message 


Error message 


I = -l 


Compute the 
function value 


Function 
value = 1 


If J is an odd 
number, function 
value = —1 
If J is an even 
number, function 
value = 1 


K-l 


Compute the 
function value 


Function 
value = 1 


Function 
value = 



Statistical results: 

f(x)-g(x) 



M(c)=Max 



f(x) 



Maximum relative error 
produced during testing 



f(x,)-g(x,) 



«x.) 



M(E)=Max | f(x)-g(x) 



a Standard deviation (root- 
mean-square) of the rela- 
tive error 

Maximum absolute error 
produced during testing 



f(x,)-g(x.) 



Standard deviation (root- 
mean-square) of the ab- 
solute error. 



In the formulas for the standard deviation, N repre- 
sents the total number of arguments in the sample; i is 
a subscript that varies from 1 to N. Appendix F ex- 
plains other symbols used above. 

Test ranges, where they do not cover the entire legal 
range of a subroutine, were selected so that users may 
infer from the accuracy figures presented the trend of 
errors as an argument moves away from the principal 
range. The accuracy of the answer deteriorates sub- 
stantially as the argument approaches the limit of the 
permitted range in several of the subroutines. This is 
particularly true for trigonometric functions. However, 
an error generated by any of these subroutines is, at 
worst, comparable in order of magnitude to the effect 
of the inherent rounding error of the argument. 



Table 4. Exponentiation With Real or Double-Precision Base 
and Integer Exponent 



Base (A) 


Exponent (J) 


J>o 


J = o 


J<o . 


A>0 


Compute 
function value 


Function 
value = 1 


Compute 
function value 


A = 


Function 
value = 


Error message 


Error message 


A<0 


Compute 
function value 


Function 
value = 1 


Compute 
function value 



Table 5. Exponentiation With Real or Double-Precision Base 
and Exponent 



Base (A) 


Exponent (B) 


B>0 


B = 


B<0 


A>0 


Compute 
function value 


Function 
value = 1 


Compute 
function value 


A = 


Function 
value = 


Error message 


Error message 


A<0 


Error message 


Function 
value = 1 


Error message 
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Section 2: Service Subprograms 



The service subprograms supplied with Fortran iv are: 

Pseudo sense light subprograms (SLITE, SL1TET) 

STOP, EXIT, and PAUSE subprograms 

Dump subprograms (DUMP, PDUMP) 

Overflow and underflow subprograms (OVERFL, DVCHK) 

Specification exception subprograms 

These subprograms are briefly described below and 



in Table 6. In most cases the actual entry point name of 
the subprogram is identical to the command name. 
However, when the user keys in the EXIT, STOP or 
PAUSE command, the compiler translates the com- 
mand name into a separate entry point name to call the 
subprogram. Both names are shown in Table 6. Further 
information concerning their usage is given in IBM 
FORTRAN IV. 



Table 6. Summary of Service Subprograms Characteristics 
1 2 



Function 


Entry Name 


Error 
Condition 


Storage 
Estimates 


Module 

Name 


HEX 


DEC 


Pseudo sense 
light sub- 
programs 


Turn all sense 
lights off or one 

sense light on 


SLITE 


Argument 
other than 
0, 1, 2, 3, 4 


324 


804 


CHCBE 


Test a sense 
light or record 

its status 


SLITET 


Argument 
other than 
1, 2, 3, 4 






CHCBE 


Overflow 
and under- 
flow sub- 
program 


Test and record 
status of 

exponent over- 
flow and under- 
flow indicators 


OVERFL 








CHCBE 


Divide cheek 
subprogram 


Test and record 
status of divide 

check indicator 


DVCHK 








CHCBE 


Exception 
processing 
subprograms 


Process arith- 
metic exceptions 


CHCBE3 ( exponent overflow ) 

CHCBE4 ( exponent underflow ) 
CHCBE5 (divide check) 








CHCBE 


Process specifi- 
cation 
exceptions 


CHCBE2 (specification) 








CHCBE 


Exit sub- 
program 


Terminate 

execution 


EXIT (CHCIW1) 
STOP (CHCIW2) 
PAUSE (CHCIW3) 




1AC 


428 


CHCIW 


Dump sub- 
program 


Dump specified 
storage area 
with or without 
termination 


DUMP, PDUMP 




48 


168 


CHCIV 



Pseudo Sense Light Subprograms 

The program-simulated machine indicator subpro- 
grams test the status of pseudo indicators, and return 
a value indicating the result of this test to the calling 
program. When the indicator is 0, it is off; when the 
indicator is other than 0, it is on. In the following de- 
scriptions of the subprograms, i represents an integer 
expression, and / represents an integer variable. 

The call slite statement is used to alter the status 
of pseudo sense lights; the call slitet statement is 
used to test, and/ or record their status. The particular 
user reference name used in the call statement de- 
pends upon the operation to be performed. 

slite is used if the four sense lights are to be turned 
off or one sense light is to be turned on. The source- 
language statement is 

CALL SLITE(t) 

where i has a value of 0, 1 ? 2, 3, or 4. 

If the value of i is 0, the four sense lights are turned 
off; if the value of i is 1, 2, 3, or 4, the corresponding 
sense light is turned on. If the value of i is not 0, 1, 2, 
3, or 4, error message 216 is issued, and execution is 
terminated. 

slitet is used if a sense light is to be tested and its 
status recorded. The source-language statement is 

CALL SLITET (f,/) 

where i has a value of 1, 2, 3, or 4, and indicates which 
sense light to test; / is set to 1 if the sense light is on or 
to 2 if the sense light is off. 

If the value of i is not 1, 2, 3, or 4, error message 216 
is issued and execution is terminated. 



DUMP and PDUMP Subprograms 

The call dump and call pdump statements allow the 
user to request that data contained within his program 
be dumped in one of nine formats. The dumps pro- 
duced will be added to the users sysout. 

It is also possible to obtain dumps using the facilities 
of the Program Control System (pes). For information 
concerning pes, see FORTRAN Programmers Guide 
and Command System User's Guide. 

The call dump statement is used if execution is to 
be terminated after the dump is taken. The source-lan- 
guage statement is 

CALL DUMP (a l9 b lt f l9 . . . , fl n ,b n ,f n ) 

where a and b are variables that indicate the limits of 
storage to be dumped ( either a or b may represent the 
upper or lower limits of storage). The clump format is 
indicated by / and may be one of the integers given in 
Table 7. A sample printout for each format is given in 
Appendix D. 



If execution of the object module is to be resumed 
after the dump is taken, the call pdump statement is 
used. The source-language statement is 

CALL PDUMP (aubufu . . . , a m b m f n ) 

where a, b, and / have the same meaning as explained 
previously. 

Table 7. DUMP/PDUMP Format Specifications 



Integer 


Specified Format 





hexadecimal 


1 
2 
3 
4 
5 


logical *1 
logical *4 
integer *2 
integer *4 
real *4 


6 


real *8 


7 
8 
9 


complex *8 
complex *16 
literal (character) 



Programming Considerations 

1. If the format control integer f is omitted, it is as- 
sumed to be equal to 0, and the dump will be hexa- 
decimal. 

2. The arguments a and b should be defined in the 
program in which the dump or pdump statement oc- 
curs; otherwise, the compiler will assign arbitrary 
addresses to them. 

3. If the program in which dump or pdump occurs is 
a subprobram, and if a and b are argument names, 
a range of storage from the calling program will be 
dumped. However, if one is an argument name and 
the other is not, unpredictable and probably large 
areas of storage will be dumped; this should be 
avoided. 

4. If one of the limits (a or b) of storage definition 
variable names is in common and the other is not 
or if it is a different (named) common, unpredict- 
able and probably large areas of storage will be 
dumped; this situation should be avoided. 

5. The literal format in Table 7 causes the area that is 
to be dumped to be treated as a string of alpha- 
meric characters. 

STOP, EXfT, and PAUSE Subprograms 

The stop, exit, and pause subprograms are called by 
the compiled object programs as a result of the source 
statements 

CALL EXIT 

STOP 

PAUSE 
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Statements that cause the users program to be termi- 
nated are 

CALL EXIT 
STOP 

If stop is issued in a conversational task, a message 
is written on the user's terminal, and control is returned 
to the terminal for entry of the next command by the 
user. If stop is issued by a nonconversational task, the 
message is written ,on the sysout data set, and the next 
command is taken from the sysin data set. The stop 
statement has the same effect when used in either a 
subprogram or main program. The call exit state- 
ment is equivalent to a stop statement. 

A pause statement executed in a program running in 
a nonconversational task will result in any associated 
messages being written to sysout; the program then 
continues execution. In a conversational task the sys- 
tem prints, at the terminal, the word pause followed by 
00000 or a l-to-5-digit integer constant, or a message, 
depending on how the operand field of the pause 
statement was written. The system then transfers con- 
trol to the terminal and awaits the user's input before 
resuming program execution. 

Overflow and Underflow Subprograms 

The call overfl statement allows a test for prior oc- 
currence of an exponent overflow or underflow excep- 



tion. The value returned by this call indicates which 
of these two conditions occurred last. After testing, the 
overflow or underflow indication is no longer available. 
The source language statement is 

CALL OVERFL (/) 

where / is set to 1 if a floating point overflow condition 
(i.e., — 16 r,:t ) exists; is set to 2 if no overflow or under- 
flow condition exists; or to 3 if a floating point under- 
flow* (i.e., <16"" <5r> ) condition exists. A more detailed 
description of each exception is given in Appendix E. 



Divide Check Subprogram 

The call dvchk statement allows a test for prior oc- 
currence of a floating point divide-check exception, 
and returns a value that indicates the existing condi- 
tion. ( Fixed-point divide checks are ignored by for- 
TRAN-compiled programs.) After testing, the indication 
of a prior divide check is no longer available. The 
source-language statement is 

CALL DVCHK (/) 

where / is set to 1 if the divide-check indicator was on, 
or to 2 if the indicator was off. A more detailed de- 
scription of the divide-check exception is^given in Ap- 
pendix E. 
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Section 3: I/O Subprograms 



This section discusses the functions, entry require- 
ments, error checks, and data references of the tss/360 
Fortran i/o library in executing the Fortran i/o state- 
ments. READ, WRITE, REWIND, BACKSPACE, END FILE, 

print, and punch. 

This section is written for both Fortran and assem- 
bler-language programmers. The Fortran programmer 
may be interested in the assumptions that the i/o rou- 
tines make, the error conditions that they check for, 
and the actions they take in case of error. The assem- 
bler-language programmer may be interested in the 
advantages of Fortran i/o facilities, particularly the 
data conversion, list-processing, and DCB-maintenance 
routines. The assembler-language programmer should 
read this section after reading IBM Time Sharing Sys- 
tem: FORTRAN Programmers Guide, Form C28-2025, 
"Appendix E. Specification of Data Set Characteris- 
tics," and IBM Time Sharing System: IBM FORTRAN 
IV, Form C28-2007, the sections titled "Input/Output 
Statements," and "Elements of the Language." Of the 
section on elements of the language, he need only 
read the subsections titled "Constants," "Variables," 
and "Arrays." 



Overview of the FORTRAN I/O Library 

There are twenty-one Fortran i/o routines. Only three 
routines, Control Initialization ( chcia ) , List Item Proc- 
essor (chcie), or List Termination (chciu), can take 
control from, or return control to, a Fortran object 
program. Thus, the Fortran i/o library can be re- 
garded as three subprograms and a number of sub- 
routines of these subprograms. 

Since the assembler-language programmer has tech- 
niques (described in Appendix B) for linking to any of 
the Fortran i/o routines, he can look upon any one of 
these routines as a subprogram. 

Another way of looking at the Fortran i/o Library 
is as two main categories of routines: i/o language 
control routines and data conversion routines. The rou- 
tines of each group interact with one another bv means 
of a common communication and work region in a 
common psect. 



I/O Language Control Routines 

There are two types of i/o language control routines: 
i/o operation control and i/o list control These rou- 
tines analyze the user's i/o requests to determine in- 
formation such as: the type of i/o operation to be per- 
formed; the number and type of list items present, if 
any; the type of format control, if any; and the i/o 
statement relationships with a user-specified ddef 
command. 

I/O Operation Control Routines 

These routines control the i/o request by creating, if 
necessary, a data control block (dcb), and analyzing 
format and namelist control specified by the user. 
After this information is processed, the i/o operation 
control routines interface with the tss data manage- 
ment routines that actually fulfill the i/o request. The 
interface with data management is accomplished by 
the routines chcib and choc, via the data manage- 
ment macro instruction facilities. 

f/O List Control Routines 

These routines examine the list items, if any, in each 
i/o request to determine the type of conversion to be 
performed. After the type of data conversion is deter- 
mined, control is given to the i/o operation control 
routines which in turn call the appropriate data con- 
version routines for final processing. 

Data Conversion Routines 

The data conversion routines are subdivided into rou- 
tines used for input processing and routines used for 
the preparation of output. These routines can process 
all the permissible types of FORTRAN-formatted data 
specified in either a format or namelist statement. 
When converting a user's data, the data conversion 
routines interact with each other according to the re- 
quirements of the user-specified format or namelist 
control. For example, for input data that is defined by 
a G-format conversion code, the General Input Con- 
version routine (ciicis) is called. This routine analyzes 
the data type to determine whether it is integer, real, 
logical, or alphameric and calls the appropriate data 
conversion routine. 
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I/O Language Control Routines 
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Figure 1. Functional Flow of FORTRAN I/O Routines 
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Routine Interrelationships 

Table 8 presents the calling relationships between the 
user program, the Fortran i/o routines, Data Manage- 
ment, and the Supervisor. 

Table 8. Calling Relationships of I/O Routines 
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The following figures describe the relationships be- 
tween routines when fulfilling a particular i/o opera- 
tion. Since the relationships vary, depending on the 
kind of i/o operation being performed, a separate dia- 
gram is presented for each of the basic i/o operations. 
Exceptions to the logical flows presented in this sub- 
section are described in detail under the individual 
routine descriptions in the following subsection. 

The type of i/o operation and its related figure ref- 
erence are; 



TYPE OF OPERATION ( FUNCTION ) 

Formatted READ with List 
Formatted READ without List 
READ with NAME LIST 
Unformatted READ with List 
Unformatted READ without List 
Formatted WRITE with List 
Formatted WRITE without List 
WRITE with NAMELIST 
Unformatted WRITE with List 
Unformatted WRITE without List 
REWIND, BACKSPACE, END FILE 



FIGURE 
2 

3 

4 

5 

6 

7 

8 

9 
10 
11 
12 
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Figure 12. BACKSPACE, REWIND, and END FILE 



Routine Descriptions 



This subsection identifies the functions, attributes, en- 
tries, routines called, error checks and data references 
of each of the twenty-one Fortran i/o routines. The 
assembler-language user should read this subsection 
in conjunction with Appendix B. 

Certain information is common to most routines; this 
information includes: a description of the attributes of 
each routine and parameter-list formats common to 
data conversion routines. 



Attributes 

Unless otherwise stated, all Fortran i/o routines are 
nonprivileged, reenterable, closed routines residing in 
syslib. ciicia, chcie, ciiciu, chciw ( except at CHCIW4), 
ciiciv, and chcbd are entered by standard Type I link- 
age with the address of a parameter list in register 1, 
and exit is a return to the calling routine. All of the 
other i/o routines are entered by restricted Type IV 
linkage. Unless otherwise stated in the description for 
a given routine, all routine exits are assumed to be 
returns to the calling routine. 



Data Conversion Routines' Parameter List 

All of the data conversion routines have a common 
parameter list in the i/o common psect. Certain data 
conversion routines do not use all the fields of the 
parameter list, in which case the fields are set to zero. 
Table 9 shows the format of the data conversion 
routines* parameter list and indicates the fields sup- 
plied by the appropriate data conversion routine. Note 
that in some cases the parameters are supplied as part 
of a common setup but either are not used by the rou- 
tine itself or are used only to pass on as parameters to 
other i/o routines. 



I/O Initialization — CHCIA 

This routine is the initial Fortran i/o Library interface 
with the user. It manages the disposition of each i/o 
request by setting information switches about for- 
matted and unformatted i/o (for use by other i/o 
routines), by allocating a buffer area for output re- 



quests, and by obtaining a logical record for input 
requests. 

Every Fortran source program i/o statement gener- 
ates a call to this routine. On this call, if there is no list, 
chcia supervises the complete execution of an i/o re- 
quest. If the i/o request is a read, write, print, or 
punch with list, ciicia simply prepares the i/o librarv 
for compiler-generated calls to List Item Processor 
( chcie ) and List Termination ( cucru ) . 



Table 9. Format of Data Conversion Routines* Parameter List 



WORD 
LOCATION 


CONTENTS 


DATA CONVERSION 
ROUTINE AFFECTED 


Word 1 


Address of list item 1 


All 


Word 2 


Byte 1 Format control 
character 2 


All 


Byte 2 Scale factor 3 


CHCII and CHCIJ, 

only 


Byte 3 Scale size 4 


CHCII and CHCIJ, 
only 


Bits 1 to 4 — Byte 
size of list item, 
Byte 4 minus one 

Bits 5 to 8 — Type 
of list item 5 


All 


W r ord 3 


Address of input or output 
buffer 


All 


Word 4 


Byte size of buffer, minus 
one 


All except CHCIM 


Word 5 


Decimal fraction width 


CHCII, CHCIJ, 
CHCIN, CHCIS, 
and CHCIT only 



1. A list item is the storage area specified by a list param- 
eter in the READ or WRITE statement. 

2. G, E 5 I, F, D, L, Z, A, H, X, T, or P (N indicates 
NAMELIST). 

3. CHCII or CHCJJ tests for an EBCDIC minus sign, 
which indicates a negative scale factor. Anything else 
indicates positive scale factor. 

4. The integer preceding the 'P\ 

5. Where 01, 02, 03, and 04 represent logical, integer, 
real, and complex, respectively. 

6. The number of decimal places to the right of the deci- 
mal point. 



Section 3: I/O Subprograms 19 



Entrtj: The entry point is chciai. The parameter list 
is variable-length and has the following format: 



Word 1 


Address of a fullword containing the user- 
specified data set reference number. 


Word 2 


Address of a control byte indicating type of 
operation. 1 


Word 3 


Address of a control byte indicating whether 
a list was present in the I/O statement and 
whether any of the following parameters in 
this list are present. 2 


Word 4 
(Optional) 


Address of a FORMAT string or NAMELIST 
table. This address is included in this param- 
eter list only if the user-requested I/O opera- 
tion had an associated FORMAT or NAME- 
LIST source statement. 3 


Word 5 

(Optional) 


Address of an error exit. This address is in- 
cluded only if the user-requested I/O opera- 
tion had the ERR operand specified in his 
source statement. 


Word 6 
( Optional ) 


Address of an end-of-file exit. This address is 
included only if the user-requested I/O oper- 
ation has the END operand specified in his 
source statement. 



NOTES ; 

1. In that control bvte, READ = 128 (XW), WRITE = 
64 (X'40"), PRINT -. 32 (X'20'), PUNCH = 16 
(X'lin, REWIND = 8 (-V08'), BACKSPACE = 4 
(X'04'), END FILE = 2 (X*02'). 

2. In that control byte, the configuration is: finrdxxx, where 
f = FORMAT 'statement, 1 = LIST parameter, n = 
NAMELIST statement, r = ERR operand, d = END 
operand, and xxx bits are alwavs set to zero. Setting any 
ot the first five bits to one indicates that the correspond- 
ing elements are present. 

3. The FORMAT string is in user-written iorm, beginning 
with the first parentheses, minus the statement number 
and the word 'FORMAT', See the CIICID routine de- 
scription for the detail* of" the NAMELIST table. 

If anv optional parameter is missing, any parameters 
following it are moved up in the list and the list is 
shortened. For example, if there is no format or 
xamklist address and no error exit address, word 4 of 
the parameter list would be the end-ot-file exit address. 

Routines Called: 

• dcb Maintenance (ciicib) 

• i/o control (cncic) 

• format Processor (ciicif) 

• namklist Processor ( CIICID ) 

• PRMPT ( CZATJl) 

• Exit (ciiciw) 

Error Checks: If the user-specified data set reference 
number is negative, an error message is issued by the 
prmpt facility, and ciiciw is entered to terminate the 
user program. 



Data References: 

• Parameter lists for the modules called by this mod- 
ule. 

• A chained list of save areas to accommodate all pos- 
sible calls to other modules. 

• A table of adcons pointing to items in the work areas 
of other modules that are to be initialized. 

• The dcb prefix (generated by dcb maintenance — 
chcib) to be set with the input parameters from this 
module. 



DCB Maintenance — CHCIB 

This routine finds or initializes the data control block 
(dcb) that contains a description of the data to be 
transmitted by a user-specified i/o operation. If an 
appropriate dcb is not found, this routine allocates the 
necessary space in the dcb table and constructs a new 
dcb, including within it information about the data to 
be transmitted that the user defined in his ddef com- 
mand. 

Entry: The entry point is chcibi. The parameter list 
in the common psect is fixed-length and has the fol- 
lowing format: 



Word 1 



Address of a fullword containing the nser-supplied 
data set reference number. 



Word 2 The data set reference number. 



ciicia stores the address of the user-supplied data set 
reference number and the data set reference number 
itself, if present, in the i/o common psect. 

Routines Called: 

• Data management routines used to search for and 
read jfcb 1 (czaeb) 

• Data management routines used to allocate storage 
for dcb construction (czcga) 

• PRMPT (CZATJl) 

• Exit (ciiciw) 

Error Checks: If the user-specified data set reference 
number exceeds 99 ? an error message is issued by the 
phmpt facilitv ( czatji) and chciw is entered to termi- 
nate 1 the* user program. 

If a discrepancy exists in the user ddef command be- 
tween permissible recfm, keylen, and dsobg values, 
an error message is issued by the prmpt facility 
(czatji) and ciiciw is called to terminate processing. 
A description of the assumption Fortran i/o makes in 
initializing associated dcbs contained in Appendix C. 

1 The Job File Control Block (JFCB) is a system control block constructed 
for each data set at DDEF time. It contains information that must be 
referred by access method routines or volume mounting routines while 
the data set is OPEN, and provides a hierarchy of pointers defining 
JOBLIB, USERLIB, and SYSLIB. 
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Data References: 

• Parameter lists for the routines called by chcib. 

• Pointers to the dcb table which consists of the dcb 
Prefix, the dcb itself, and additionally two decbs if 
the user has specified Basic Sequential Access 
Method ( bsam ) in his ddef command. 

• A chained list of save areas to accommodate all pos- 
sible calls to other routines. 

Format and Content of the DCB Prefix: The dcb 
prefix is used by the Fortran i/o routines, in conjunc- . 
tion with the dcb, when performing any type of i/o 
operation. The dcb prefix, created by chcib, is eight 
words long and always immediately precedes the dcb 
itself. 

Table 10. Format and Content of DCB Prefix 



Word 1 


The address of the starting location in 
the buffer area for the current logical 
record. 


Word 2 


The address of the current location in 
the buffer area for the current logical 
record. 


Word 3 


The address of the end location in the 
buffer area for the current logical rec- 
ord. 


Word 4 


Byte 1 : 
Byte 2: 

Byte 3; 
Byte 4: 


Current operation (READ, WRITE, 

etc.) 1 

Control flags (FORMAT, NAMELIST, 

List, ERR exit, END exit) 2 

Control flags (Span, GATE, recent 
READ, END or ERR encountered) 3 
Previous operation (byte 1 from last 
call on CHCIC with this DCB ) 


Word 5 


The address of current DECS, if re- 
quired (BSAM) 


Word 6 


The user-specified data set reference 
number, plus one. 


Word 7 


The address of the next DCB. 


| Word 8 


Save area for the address of the previ- 
ous DCB for that data set reference 

number. 


i Word 9 


DCB begins here. 



notes : 

1. See parameter list at entry to CHCIA, Note 1. 

2. See parameter list at entry to CHCIA, Note 2. 

3. The configuration is; gxdrxxln, where g = GATE I/O, 
d — end of data set (END), r = error (ERR), 1 = span 
from last record, or recent READ, and n = span to next 
record. The x bits are always set to zero. All bits set to 
zero signifies that there is no span. (Spanning is used in 
the case of unformatted records, where a physical block 
size was defined. It is the process of jumping from the 
end of one record to the beginning of the next. ) 



I O Control — CHCIC 

This routine fulfills i/o requests made through other 
i/o library routines by using the data management 
macro instruction facilities of tss. The particular data 
management facilities to be used are determined 



both by the type of i/o statement issued in the user 
program, and by any related ddef commands, if any, 
defining such things as the type of records being trans- 
ferred and the manner in which they should be proc- 
essed. 

The following list identifies the more significant 
macro instructions used by chcic for each of the For- 
tran i/o statements. 



FORTAN I/O 

STATEMENT 
ftEAD 



WRITE 



REWIND 



BACKSPACE 



END FILE 



CHCIC FUNCTION 

Obtains a logical record from a user-specified 
input source by using the READ, GATRD, or 
GET macro instruction. 

Initializes the writing of a logical record by 
establishing pointers to the output buffer area. 
Subsequent output processing is performed bv 
using the WRITE, GATWR, or PUT macro 
instruction. 

Repositions the user-specified volume of one or 
more data sets to the first record of the first 
data set by using the POINT or SETL macro 
instruction. 

Repositions the user-specified data set to the 
previous logical record by using the NOTE, 
POINT, SETL, and BSP macro instructions. 
Defines the end of the user-specified data set 
by using the WRITE and STOW macro in- 
structions. 



Entry: The entry point is chcici. 

Routines Called: 

• dcb Maintenance ( chcib ) 

• Exit (chciw) 

• Data management routines to perform i/o functions 
as determined by the macro instruction issued. 

• Error message control ( chcix ) 

• FRMPT (CZATJI) 

Error Checks: If the i/o operations performed by 
data management cause either a synad 1 or eodad 2 exit, 
and if the user provided an err or end return point, 
chcic locates the adcons for these return points in the 
work area chcrww and locates the register save area 
for the user's program registers. Return is then made to 
the err or end return point rather than to the calling 
i/o routine. 

If the user did not provide return points (or if the 
operation was other than a read statement), an error 
message is issued and the program is terminated. 

If an invalid character is encountered in hexadecimal 
input from a gate 3 read operation performed for an 
unformatted read statement, an error message is issued 
and the erroneous character is treated as the termina- 
tion of the hexadecimal input. Processing then con- 
tinues. 

2 syNAD: synchronous error exit address, for automatically transferring 
control to a user-supplied routine if an uncorrectable I/O error occurs. 

2 EODAD: end of data set address, for automatically transferring control 
to an end-of-data routine when end of an input data set is detected 
during processing. 

3 GATE I/O is input from SYSIN or output to SYSOUT. 
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In addition to the above error checks, error messages 
are issued ( prmpt macro instruction ) and the user pro- 
gram is terminated by chow for any of the following 
reasons; 

• The record is not format- V for unformatted read 
statement. 

* Error return code received from the use of the find 
or stow macro instruction for a member in a vpam 
data set. 

# Invalid sequence of i/o operations for a user-speci- 
fied data set reference number. The invalid se- 
quences are: read preceded by end file; end file 
preceded by read; and read preceded bv write ( ex- 
cept when using gate i/o ) . 

Data References: 

# References to the standard dcb and its associated dcb 
prefix. 

• A chained list of save areas to accommodate all pos- 
sible calls to other routines needed. 



NAMELIST Processor— CHCID 

This routine interacts with chcic to control the i/o for 
each namelist record and interacts with the appropri- 
ate- data conversion routines to bring about the desired 
item-bv-i torn conversion. 

Kntnj: The entry point is chcidi. The parameter list 
consists of a single word; 



W< 



Address of the NAMELIST fable generated by 
the FORTRAN compiler as part of the user object 

program. 



Routines Called: 

• i/o Control (ciicic) 

• Complex Input Conversion (ciic.im) 

• Complex Output Conversion (ciicin) 

• Gen era 1 1 n pu t Con vers i on ( ci i cis ) 

• General Output Conversion (ciicit) 

• prmpt (cza'tji) 

• Exit (chctw) 

Error Checks: There are no error checks for output. 
For input, if errors are detected in the namelist table, 
a message is issued via prmpt and ciiciw is called to 
terminate the user program. Other error messages are 
generated for any of the conditions listed below. In 
these cases, processing continues with the next entry 
of the input record. 

• Name exceeds six characters 

• First character of each input record is not blank 

• Subscripts appear on a name that is not an array 
name 

• Incorrect number or range of subscripts 



• Subscripting causes array size to be exceeded 

• Multiple constants or repeated constants appear with 
a name that is not a subscripted array name, or ex- 
ceed the size of an array 

• An equal sign or left parenthesis is not preceded by 
the variable or array name for that item. 

• An invalid character appears in a repeat constant 

• End of a logical record caused an item to be logically 
incomplete 

• The namelist name is not in the namelist table. 

Data References: 

• Parameter lists for other i/o library routines called 
by this routine. 

• A chained list of save areas to accommodate all pos- 
sible calls to other routines needed. 

NAMELIST Table: The address of the namelist 
table generated by the Fortran compiler or by the 
assembler-language programmer is communicated in 
the call to i/o Initialization (chcia) and then passed 
to this routine. The table is made up of two-word en- 
tries, each of which contains an identifier in the first 
halfword. 

NAMELIST Name Entry: 
Bytes 0-1: Identifier (X'0100') 

2-7; Name (left-justified) 

Variable Name Entry; 

Bytes 0-1; Identifier (X'0200') 
2-7; Name (left- justified) 

Variable Type and Location Entry: 
Bytes 0-1: Identifier (X'0300') 

2: Length and Type ( 4 bits each ) 
Length: Number of bytes minus 1 
Type: X'01' Logical 
X'02' Integer 
X'03' Real 
X'04' Complex 
3; Class; Letter A for array; otherwise, 
an S 
4-7: Storage Location 

Array Size Entry: 

Bytes 0-1: Identifier (x'0400') 
2-3: Not used 
4-7: Number of bytes in array 

Dimension Product Entry: 
Bytes 0-1: Identifier (x'0500') 
2-3: Not used 

4-7; Dimension Product (see explanation 
below ) 

Terminal Entry: 
Bytes 0-3: Zero 

4-7: Not used 
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A dimension is a level of subdivision, or level of sub- 
scripting, within an array. For example, an array could 
be a string of seven thirty-word elements ( first dimen- 
sion), each subdivided into six five- word elements 
(second dimension), each subdivided into five one- 
word elements ( third dimension ) . An array may have 
as many as seven dimensions. 

For each dimension there is a corresponding dimen- 
sion product, which is the product of 1) the byte-size 
of the array's smallest element, 2) the number of ele- 
ments within all lower dimensions except the first di- 
mension, and 3) the number of elements within that 
dimension. In the example just given, the dimension 
product for the third dimension would be 4 x 6 x 5, or 
120. This dimension product would be seven times 
greater if there were another dimension before the 
seven-element dimension. The dimension product for 
the first dimension is always the byte-size of the array's 
smallest element — this dimension product is never en- 
tered. If there is only one level of subdivision, there 
should be no Dimension Product Entry. 

Following is a hexademical representation of the 
namelist table for a three-dimension array such as 
that described above, where the array is named 'C and 
contains real numbers. The namelist name is list. 



01 

E2 


00 

E3 


D3 

40 


C9 
40 


NAMELIST name 


02 

40 


00 

40 


C3 
40 


40 
40 


Array name 


03 
00 


00 
0E 


33 
63 


CI 
74 


Variable type 


04 
00 


00 
00 


00 
03 


00 
48 


Array size 


05 
00 


00 
00 


00 
00 


00 
18 


Dimension product 


05 
00 


00 
00 


00 
00 


00 

78 


Dimension product 


00 


00 


00 


00 


Terminal entry 



ist Item Processor — CHCIE 

Evrry i/o statement in the users source program gen- 
erates one or more calls to this routine if there is a list 
associated with a read, write, print, or punch. A list 
item may be a simple variable, an array element (a 
subscripted variable), or an entire array. If a format 
statement is specified, this routine calls on Format 
Processor (ciicif) to control any necessary conversion. 
If there is no format statement, chcie is directly re- 
sponsible for filling or emptying the output or input 
buffer area. 

Entry: The entry point is chciei. Register contains 
either zeros, if the list item is a single element, or a 
number expressing the array length, in bytes, if the list 
item is an entire array. The parameter list is fixed- 
length and has the following format: 



Word 1 


Address of a control byte. The first four bits of 




the control byte contain the size of the element, 




minus one. The second four bits contain a flag 




indicating the type of item as follows: 




Flag Type of Item 




01 logical 




02 integer 




03 real 




04 complex 


Word 2 


Address of a first (or only) element of the list 




item. 



Routines Called: 

• Format Processor (chcif) 

• i/o Control (chcic) 

• prmpt (czatjt) 

• Exit (chciw) 

Error Check: With unformatted input, if a list item is 
.requested after the logical record is exhausted, an error 
message is transmitted to the user via prmpt, and 
chciw is called to terminate the user-program. 

Data References: 

• Parameter lists for other i/o library routines called 
by chcie. 

• A chained list of save areas to accommodate all pos- 
sible calls to other routines needed. 

• A full word, CHCIB2, which is in the ciicir work area 
and contains the address of the dcb prefix. 

• The first fifteen bytes of the dcb prefix. 



FORMAT Processor— CHCIF 

This routine interacts with chcic to control the i/o for 
each FORMAT-referenced record, and interacts with the 
appropriate data conversion routines to bring about the 
item-by-item conversion specified by the format state- 
ment. 

Entry: Before the first entry to chcif to process a 
reference to a format statement, chcia ( or the assem- 
bler-language programmer, if he is bypassing chcia) 
does the following : 

• Store the address of the format character string in 
chcrww. The statement number and the word 
'format are omitted from the string. 

• Set to zero the second and third words of chcifw. 
The entry point is chcifi. The parameter list is fixed- 
length and has the following format: 



Word 1 


Address of the list item, if any. 

Zero indicates that no list item was specified. 


Word 2 


Byte size of list item and type in low order byte 
of word. (See word 1 of CHCIE parameter list.) 


Word 3 


Address of the start of the format string. 
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Routines Called: 

• i/o Control (ciicic) 

• Error Message Control (cncrx) 
® Exit (cuciw) 

• One of the eleven data conversion routines (chcih 
through ohctt) 

Error Checks: Since form at statements mav be dy- 
namically modified, certain error conditions may arise 
due to the syntax of the format string. If there are no 
syntax errors, errors could arise due to conversion of 
the data. In such cases the conversion routines issue 
messages describing the errors before returning. All 
syntax error cheeks produce messages describing the 
error. 

Processing is terminated upon encountering invalid 
control characters in the string, strings that exceed the 
maximum, or too many levels of parentheses. When it 
is possible to assume values other than those specified 
( as in the ease of invalid size of w or d fields after a 
control character), processing will continue on the cur- 
rent item after the error message is issued. Otherwise, 
the erroneous format item is skipped and processing 
continues with the next control character, 

Data References: 

• Parameter lists for the routines called by chcif. 

• A chained list of save areas to accommodate all pos- 
sible calls to other routines. 

• Counters for any repetition and scale factors encoun- 
tered. 

Integer Output Conversion— CHCIH 
Tliis routine converts a two-byte or four-byte binary 
list item to an integer field in the output buffer, accord- 
ing to the format In, where n is the integer field size. 

Entry: The entry point is ciicihj. The parameter list 
is described at the beginning of this subsection, under 
"Data Conversion Ron line Parameter Lists," 

Routines Called: 

• Error Message Control (ciicix) 

Error Checks: If the output buffer area is too small 
to contain the integer field, the field is filled with 
asterisks and a message is issued by chcix. 

Da fa Refe ren ces : 

• A parameter list for cncrx. 

• A save area to accommodate the call to ceicix. 

• A work area, chcih w, to be used by this routine. 

Reel and Integer Input Conversion— -CHOI 
This routine converts a data field in an input buffer to 
the* appropriate type list item. An integer field in the 
input buffer is converted to a binary list item. A real 
field in the input buffer is converted to a single- or 
double-precision floating-point list item. The integer 
field has a format In, where n is the field width. The 



real field has a format Fw.d, Ete.d, or Dw.d, where w 
is the field width and d is the width of the decimal 
fraction. 

Entry: There are three entry points; chcii, ciicik, 

and chcig. The parameter list is described at the be- 
ginning of this subsection under "Data Conversion 
Routine Parameter Lists." 
Routines Called: 

• Error Message Control (chcix) 

Error Checks: If the format specification (F, E, D, 
or I) is improperly specified or the data field is greater 
than the permissible range, chcix is called. 

Data References: 

• A parameter list for chcix. 

• Adcons for the table in Real Output Conversion 
( chcij ) that contains powers of ten. 

• A work area, chciiw, containing: two doubleword 
areas for calls to ciicik and chcig, and a 32-byte area 
for temporary storage. 

Real Output Conversion — CHCIJ 

This routine converts a single- or double-precision 
floating point list item to a real field in the output 
buffer. The real field has a format of either Fw.d, 
Dw.d, or Fw.d where w is the field width and d is the 
size of the fractional position, in digit positions. 

Entry: There are two entry points; chciji and 
chcili. The parameter list is at the beginning of this 
subsection under "Data Conversion Routine Parameter 
Lists." 

Routines Called: There are no calls that can occur 
besides the final return to the calling routine. 

Error Checks: If the output buffer area is too small 
to contain the real field, the real field is filled with 
asterisks. 

Data- References: 

• A table of power of ten in double-precision floating- 
point. It has an external name cnciL2, so that it can 
be referred to and used by other i/o library routines. 
The table structure is : 

CHCIL2 DC D1E1,1E2,1E3,1E4,1E5,1E6,1E7,1E8,1E9,1E10' 
DC D'lEll,lE!2aE13JEM,lE15,lE16,lE17,lE18, 
1E19.1E20' 



DC D i lE71 3 lE72 3 lE73 J lE74JE75,lE-76 > lE-77, 
1E-78* 



Complex Input Conversion — CHCIM 

This routine converts a complex data field from an 
input buffer to a complex list item, consisting of two 
real data fields. Each real field is converted to a single- 
or double-precision floating-point list item according 
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to the format Fw.d, Fw.d, or Dw.d, where w is the 
real field width and d is the width of the decimal 

fraction. 

Entry: The entry point is chcimi. The parameter list 
is described at the beginning of this subsection under 
"Data Conversion Routine Parameter Lists." 

Routines Called: 

• Real and Integer Input Conversion (chcii) 

• PRMPT (CZATJI) 

Error Checks: If only one or if no real fields exist in 
the complex data field in the input buffer, or if there is 
a missing parentheses or central comma, chcim issues 
an error message via prmpt. No further action is taken 
and the list items remain unchanged. If either or both 
real fields contain invalid characters or exceed the per- 
missible magnitude range, chcii assumes the responsi- 
bility for producing an error message. 

Data References: 

• Parameter lists for routines called by chcim. 

• Adcons for the table produced by chcij, containing 
powers of ten. 

• A chained list of save areas to accommodate all pos- 
sible calls to other routines. 



Complex Output Conversion — CHCIN 

This routine converts a complex list item consisting of 
two, single- or double-precision floating point items to 
a complex data field in an output buffer. Each floating 
point list itern is converted to a real data field accord- 
ing to the format code Fw.d, Fw.d, Dw.d, or Giv.s, 
where to is the real field width, d is the width of the 
decimal fraction, and s is the number of significant 
digits. 

Entry: The entry point is chcixi. The parameter list 
is described at the beginning of this subsection under 
"Data Conversion Routine Parameter Lists." 

Routines Called: 

• General Output Conversion (chcit) 

Error Check: If the format specifications (F, E ,D, 

or G) is improperly specified or the real data field is 
greater than the permissible range, the general output 
conversion routine (chcit) assumes the responsibility 
for producing an error message. 

Da ta Refe re nces: 

• P:'.: ameter list for chcit. 

• Adcons for the table produced by chcij, containing 
powers of ten. 

• A chained list of save areas to acommodate all pos- 
sible calls to other routines needed. 



Alphameric and Hexadecimal Input 
Conversion — CHCIO 

This routine transfers a specified number of bytes 
( alphameric or hexadecimal characters ) from an input 
buffer area to a list item. The format is Aw (alpha- 
meric) or Zw (hexadecimal), where w, field width, is 
the number of characters being transferred. 

Entry: The entry points are chcioi (alphameric data) 
and CHC102 ( hexadecimal data ) . The parameter list is 
described at the beginning of this subsection under 
"Data Conversion Routine Parameter Lists." 

Routines Called: None. 

Error Checks: None. 

Data References: None. 

Alphameric and Hexadecimal Output 
Conversion — CHCIP 

This routine transfers a specified number of bytes ( al- 
phameric or hexadecimal characters) to an output 
buffer area from a list item. The format is Aio (alpha- 
meric) or Zw (hexadecimal), where w, field width, is 
the number of characters being transferred. 

Entry: The entry points are chcipi (alphameric 
data) and chcip2 (hexadecimal data). The parameter 
list is described at the beginning of this subsection 
under "Data Conversion Routine Parameter Lists." 

Routines Called: None. 

Error Checks: None. 

Data References: None. 

Logical input Conversion — CHCIQ 

This routine converts a logical field in the input buffer 
area. The logical field has the format hie, where to is 
the logical field width. 

Entry: The entry point is chciqi. The parameter list 
is described at the beginning of this subsection under 
"Data Conversion Routine Parameter Lists." 

Routifies Called: None. 

Error Checks: None. 

Data References: None. 

Logical Output Conversion — CHCIR 

This routine converts a list item to a logical field in the 
output buffer area. The logical field has the format Lw, 
where w is the logical field width. 

Entry: The entry point is chciri. The parameter list 
is described at the beginning of this subsection, under 
"Data Conversion Routine Parameter Lists." 

Routines Called: None. 

Error Checks: None. 

Data References: None. 

General Input Conversion — CHCIS 

This routine converts a data field in the input buffer to 
a list item according to the format Gw.s, where w is 
the field width and s is an optional specification of the 
number of significant digits. 
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Entry; The entry point is chcisi. The parameter list 
is described at the beginning of this subsection under 
"Data Conversion Routine Parameter Lists." 

Routines Called; 

• Real and Integer Input Conversion (chcii) 

• Logical Input Conversion ( chciq ) 

• Alphameric Input Conversion (ciicio) 

Error Checks; chcis performs no error checking. 
Error checks, if any, are made by the called data con- 
version routines. 

Data References; 

• Parameter lists for the routines called by chcis. 

• A chained list of save areas to accommodate all pos- 
sible calls to other routines. 



General Output Conversion — CHCIT 

The routine converts a list item to a data field in the 
output buffer, according to the format Gw.s, where w 
is the field width and s is an optional specification of 
the number of significant digits. 

Entry; The entry point is chciti. The parameter list 
is described at the beginning of this subsection under 
"Data Conversion Routine Parameter Lists." 

Routines Called; 

• Integer Output Conversion (chcih) 

• Real Output Conversion (chcij) 

• Logical Output Conversion (ciicir) 

Error Checks: chcit performs no error checks. Dis- 
crepancies between the size and type specification of 
the list item and the data field are detected by the 
called conversion routine. 

Data References: 

• Parameter lists for the routines called by chcit. 

• A chained list of save areas to accommodate all pos- 
sible calls to other routines. 



List Termination — CHCiU 

This routine terminates list processing for a read, 
write, print, or punch statement, and completes any 
i/o operation that is pending. 

Entry: The single entry point is ciiciui. No param- 
eters are passed. 

Routines Called; 

• Format Processor (chcif) 

• i/o Control (chcic) 

The final return is made with registers unchanged, 
except that register 13 will be set to the address of the 
calling module's psect and register 15 will be set to 
zero. 

Error Check: None 

Data References; 

• Parameter lists for other i/o library routines called 
by chciu. 



• A chained list of save areas to accommodate all pos- 
sible calls to other routines. 

• A control byte within the dcb prefix that describes 
the current operation. (See "DCB Maintenance — 
CHCIB" and Table 10, "Format and Content of the 
DCB Prefix," in this section. ) 



Exit— CHCIW 

The Exit Routine's subprograms, stop, exit, and pause, 
are described in Section 2. 



Error Message Control — CHC1X 

This routine receives the text of error messages from 
other i/o library routines during execution, and deliv- 
ers those messages as output via the gate macro in- 
struction, to the user's sysout. In conversational mode, 
for example, any error message generated is passed to 
this routine for transmission to the user's terminal. 

Entry: The entry point is chcixi. The parameter list 
is fixed-length and has the following format; 



Word 1 


Address of first part of message. 


Word 2 


Byte length of first part of message, minus one. 


Word 3 


Address of second part of message. 


Word 4 


Byte length of second part of message, minus one. 



The first part of each message is a character string 
that never changes for that message, and is therefore 
part of the calling routine's csect. The second part is 
some data item that does change ( such as the contents 
of a field containing invalid characters), and which, 
therefore, is in a psect (either of the user's problem 
program or of the i/o library routines). If only a single 
part message is to be transmitted to sysout, word 3 of 
the parameter list is set to zero. 

Routiries Called: 

• gatwr macro instruction 

Error Checks; The size of the second part of a mes- 
sage must not be greater than 49 bytes. If this limit is 
exceeded, only the leftmost 49 bytes of data will be 
obtained from the invalid field. No error message is 
generated for this situation. 

Data References; 

• A 100-byte buffer area used for the error message. 

• Parameter lists for the routines called by chcix. 

• A chained list of save areas to accommodate all pos- 
sible calls to other routines. 
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Interruption and Machine Indicator Routine — CHCBD 

This routine sets bits in the psw so that the fixed- 
point overflow and significance exceptions will be ig- 
nored, and directs the system interruption handler 
where to pass control if any of the following four 
exceptions occur: 

Exception Subprogram 

Specification chcbe2 

Exponent overflow chcbe3 

Exponent underflow CHCBE4 

Divide check chcbes 

In addition, this routine initializes the machine indi- 
cator flags and the sense light indicators, and clears 
any pointers to entries in the dcb table. It then returns 
control to the calling program. 

Entry: The entry point is chcbdi. There are no entry 
parameters. 



I/O Communication — CHCIY 

This table contains space for linking register save 
areas and an area in which to construct a chain of 
DCBs. 

The format of chciy within the i/o psect communi- 
cation region (i.e., save and dcb areas) is; 



PSECT Communication Region 



76 



152 



228 



304 



380 



456 



SAVE 1 


SAVE 2 


SAVE 3 




SAVE 4 


SAVE 5 


SAVE 6 


460 464 468 


472 


CHCIY9 



6 19-Word Save Areas 
( Each area has the 
address of the next 
save area in word 19 ) 



Pointers to work areas 
for CHCIA, CHCIE, 
CHCIU, CHCIV, 
CHCIW 

Area for construction 
of DCBs and DCB 
prefixes 
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Appendix A: Replacing FORTRAN IV Library Subprograms 



This appendix provides a general description of tech- 
niques for replacing a Fortran iv library subprogram 
with a "private" version of the same program. The dis- 
cussion below does not describe a technique for replac- 
ing the copy of a subprogram in a manner that will 
cause all users of Fortran iv library subprograms to 
use the new version. 

It is recommended that a user-written version be 
loaded explicitly, with a load command. The Fortran 
iv mathematical subprograms, service subprograms, 
and i/o subprograms reside in syslib as six link-edited 
modules, and implicit loading of a user-written version 
is possible only when the corresponding Fortran iv 
library module is not already loaded. 

Many subprograms call other subprograms, as shown 
in Table 1, Table 8, and Figures 1-11. For example, the 
csqrt subprogram, called by a Fortran program to 
find the square root of a complexes number, requires 
the cabs and sqrt subprograms. If the Fortran user 
loaded his own version of sqrt, the csqrt subprogram 
would use this version. Note that if the Fortran user 
wishes the csqrt subprogram to use his own version of 
sqrt, he must supply the entire mathlib (since it is 
link-edited). The user may not supply one routine only 
without performing a new link-edit. 

The Fortran compiler and the Fortran iv library 
subprograms expect a substituted subprogram to sat- 



isfy the same references as the original subprogram. 
The following table serves as a guide to the external 
names of each subprogram. 

Table 1 1 . External Names of FORTRAN IV Library 
Subprograms 





MATHEMATICAL 


SERVICE 


I/O 




SUBPROGRAMS 


SUU PROGRAMS 


SUBPROGRAMS 


Entry 


See Tables 1 


See Table 6. 


See Section 3. 


Name 


and 2. 






Routine 


See Tables 1 


See Table 6. 


See Section 3. 


Name 


and 2. 






CSECT 


Routine name 


CHCBD and 


Routine name 


Name 


suffixed by *W\ 


CHCBE: Rou- 


suffixed bv *C\ 






tine name 


CHCIR and 






suffixed by 'W\ 


CHCIC have 






CHCIV and 


additional 






CHCIW: Rou- 


CSECTS with 






tine name 


routine name 






suffixed by *C\ 


suffixed by *X\ 


PSECT 


Routine name 


CHCBD and 


Routine name 


Name 


suffixed by *R\ 


CHCBE: Rou- 
tine name 
suffixed bv *R\ 
CHCIV and 
CHCIW: Rou- 
tine name 
suffixed by *W\ 


suffixed by 'W. 
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Appendix B: Assembler Language Information 



The mathematical, service, and i/o subprograms are 
available to the tss assembler-language programmer. 
The following explains the method of calling a library 
subprogram from an assembler-language program and 
gives other information for the assembler-language 
programmer who wants to use these subprograms. 
Before reading any subdivision of this appendix, the 
assembler-language programmer should become fa- 
miliar with the corresponding section of the main text. 

Note: The examples in this appendix have not been 
tested on the current system. 

The linkage from Fortran compiled programs to 
Fortran rv subprograms is a standard, Type I linkage. 
Assembler-language programmers must link to these 
subprograms using an identical linkage. The call 
macro instruction provides a number of different means 
for establishing the correct linkage. (See Assembler 
User Macro Instructions.) fy hand-coded linkage may 
also be used, but such linkages should generally be 
avoided when macro instructions supply the service re- 
quired. Regardless of which form of linkage is used, 
however, the register usages for linkage are: 

1. Register 1 must point to whatever parameter list the 
subprogram requires. 

2. Register 13 must point to a 19-word save area in the 
calling program. 

3. Register 14 must contain the address in the calling 
program to which control will be returned by the 
called program at the completion of its operation. 

4. Register 15 must be loaded with the address of the 
entry name, and this register is used to transfer con- 
trol to the called program. 

Before returning to the calling program, Fortran 
librarv subprograms always restore general registers 1 
through 14. General register is restored except 
when the result is returned by a mathematical sub- 
program and is an integer, in which case the integer is 
contained in this register. The floating registers are not 
restored, and should be assumed destroyed. General 
register 15 is not restored, as future modifications to 
the Fortran library subprograms may make use of this 
register for a return code (they do not currently do 
so ) ; this register should be assumed destroyed. 

Mathematical Subprograms 

The parameter list for a mathematical subprogram 
must contain the addresses of the arguments in the 
proper order: 



• Directly referenced subprograms. The order is the 
same as that in the list of operands within the paren- 
theses in the corresponding Fortran source state- 
ment. For example the source statement 

ANS = SIN (RADIAN) 

in Fortran coding corresponds to an assembler-lan- 
guage call containing one address in the parameter 
list — the address of radian. The Fortran statement 

ANS = ATAN2(X,Y) 

produces a linkage with a parameter list containing 
the addresses of X and Y, in that order. The assem- 
bler-language programmer's linkage to atan2 must 
do the same. 

• Indirectly referenced subprograms. The order for the 
exponentiation subprograms is; address of the num- 
ber to be raised to a power and the address of the 
power itself. 

The arguments pointed to by the parameter list can 
be either integer values, or normalized floating point 
real or complex values, as required by the called pro- 
gram. An integer argument occupies four locations of 
storage. A real argument occupies either four or eight 
locations of storage. An argument occupying eight lo- 
cations of storage starts on a doubleword boundary 
and occupies two adjacent words. The address of the 
first word is the address of the entire argument. 

A complex argument occupies either eight or six- 
teen locations of storage, starts on a doubleword bound- 
ary, and occupies adjacent words. The first half of the 
argument contains the real part of the complex argu- 
ment; the second half contains the imaginary part. The 
address of the real part of the argument is the address 
of the entire argument. 

Each mathematical subprogram returns a single 
answer — either an integer value, a normalized floating 
point value, or a complex value. An integer answer is 
stored in general register 0, a real answer is stored in 
floating point register 0, and a complex answer is stored 
in floating point registers and 2. The real and com- 
plex parts of a complex number occupying eight stor- 
age locations will be in the high-order four storage 
locations of floating point registers and 2. 

Examples of the use of the CALL macro instruction 
for an assembler-language programmer using the sine 
program are: 
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LA 

CALL 

LA 
LA 
CALL 



SAVE DS 

PARLIST DC 
VCON ADCON 



13, SAVE 

SIN, (RADIAN) 

13, SAVE 

15, VCON 

(15). MF=(E, PARLIST) 



19F 

A(RADIAN) 

IMPLICIT, EP=SIN 



Point to a 19- 
word save area. 



The above examples produce code equivalent to the 
following hand-coded linkages. ( Several additional in- 
structions are included for greater clarity. ) This exam- 
ple assumes that appropriate cover registers have been 
established, and radian - contains the value, in radians, 
for which the sine is to be obtained. 

Point to a 19-word save area. 

Point to the parameter list. 

Store the R-con in the 19th 
vvorcl of the callers save area. 

Obtain the address of the entry 
point. 

Branch to the entry point, set- 
ting register 14 to the address 
c5f the instruction following 
the BASR. 
STE 0, ANS Store the result in ANS. 



LA 


13, 


SAVE 


LA 


1, PARLIST 


L 


14, 


RCON 


ST 


14, 


72 (a 13) 


L 


15, 


VCON 


BASR 14, 


15 



SAVE 
PARLIST 

VCON 
RCON 

RADIAN 

ANS 



DS 
DC 

DC 
DC 

DS 
DS 



19F The 19-word save area. 

A (RADIAN) The sine at RADIAN is to be 

computed. 
V (SIN) The V-R-con pair for the 

R (SIN) system entry to the sine 

program. 



F 
F 



The result is stored here. 



Service Subprograms 

The calling sequence for dump and pdump may specify 
a variable number of parameters. Forms of the call 
macro instruction are available for this purpose. The 
linkage is identical to that described above, with one 
exception: immediately preceding the address of the 
first parameter there must be a word containing, in 
binary and right adjusted, the number of addresses in 
the parameter list. Note that this word contains a 
count, not the address of a count. 



I/O Subprograms 

As with other i/o, data sets used with the Fortran 
i/o library must be defined. Unless the program is 
using gate i/o, the programmer must give a ddef com- 
mand. For example: 

DDEF DDNAME=FT10F001,DSORG=VS,DSNAME=PAY 

This command is presented in keyword form, for clar- 



ity. It could also be written in the shorter, positional 
form as follows : 

DDEF FT10F001,VS,PAY 

Note that the ddname is in Fortran format and con- 
tains the data set reference number in the two digits 
following the 'ft/ 

Having satisfied ddef requirements, the programmer 
is in position to implement the information given in 
Section 3: I/O Subprograms, The following arc exam- 
ples of ways the assembler-language programmer 
might use Fortran i/o facilities. 

Formatted READ with List 

Assume that the programmer wants to read an eighty- 
byte record containing three integer numbers in the 
first half of the record. The first number occupies bytes 
three through eight, the second occupies bytes fifteen 
and sixteen, and the third occupies bytes thirty-nine 
and forty. The rest of the first forty bytes are blank. 
The second forty bytes are to be ignored. 

The numbers are to be converted from character to 
integer form and placed in storage areas (list items) 
labeled A, B, and C, respectively. 

The programmer chooses not to construct a dcb, 
since chcib (dcb Maintenance) will construct one for 
him when it finds that there is no dcb for the data set 
reference number given in the ddef command. 



LA 
CALL 



13,SAVE 

CHCIA1, (PARLISTO) 



The linkage shown by arrow 
4 of Figure 2, to CHCIA 
( I/O Initialization ) . 

At this point, chcia ( 1 ) causes chcib to create the dcb, 

(2) causes chcic (i/o Control) to perform the i/o, and 

(3) passes the format string to chcie (List Item 
Processor ) . 

SR 0,0 Indicate to CHCIE that the 

list item is a single 
element. 

CALL CHCIE1, ( PARLIST1 ) The linkage shown by arrow 

2 of Figure 2, to CHCIE. 
CHCIE will process the 
first list item. 



The second list item. 



SR 


o ? o 


CALL 


CHCIE1, (PARLIST2) 


SR 


0,0 


CALL 


CHCIE1,(PARLIST3) 


CALL 


CHCIU1 



The third list item. 
The linkage shown by arrow 
7 of Figure 2, to CHCIU 
( List Termination ) . 
There are no parameters. 



SAVE DS 19F 

* PARAMETER LIST FOR CHCIA 
PARLISTO DC A(DSRN) 
DC A(CREAD) 
DC A(COPNDS) 
DC A (FORMAT) 
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DC 
DC 



A(LABELl) 
A(LABEL2) 



CREAD DC 



COPNDS DC 



Addresses of the user- 
written error-handling 
and end-of-file routines. 
Both parameters are 
optional. 
DSRN DC XL4'0A* The data set reference 

number (10io). The I/O 
routines expect it to be in 
fullword, binary form. 

X'80 The control byte addressed 

by the second word of 
the parameter list. 
Signfiies READ 
operation. 

X*D8' The control byte addressed 

by the third word of the 
parameter list. Specifies 
that there will be list 
processing, and that there 
are entries in the last 
three words of the 
parameter list. 
DS OF Puts FORMAT string on a 

fullword boundary. 

Following is the format string. Note that the fields 
are defined in such a way that the numbers are in the 
rightmost portions of the fields. This must always be 
done with integer conversion, since blanks are treated 
as zero and would multiply any integer value by ten 
for every blank on the right. 

The FORMAT string. 
CHCIE 



FORMAT DC 


C'(G8,G8, 




G24)* 


* PARAMETER 


LISTS FOR 


PARLIST1 DC 


A (ITEM) 


DC 


A(A) 


PARLIST2 DC 


A (ITEM) 


DC 


A(B) 


PARLIST3 DC 


A(ITEM) 


DC 


A(C) 


T rT \: DC 


X'32' 



The first four bits of this 
control byte indicate that 
the list item ( into which 
an integer will be placed) 
is four (3+1) bytes long. 
The second four bits 
indicate that the 
characters which the 
FORMAT statement 
causes to be read are to 
be converted into integer 
form. 



FORMAT Conversion and Lis* Processing 
Assume that the progra.ir.mer has scanned numbers 
into hold, a 400-byte area. The numbers are in ebcdic 
form, with the format xxx.xxx, where x" is any digit. 
They occupy contiguous, two-word elements. The pro- 
grammer wants to convert them into real form and 
move the result into a 50-word array. (An array is 
simply a string of equal-length elements.) The pro- 
grammer wants to use the Fortran i/o library only for 
its data conversion and list processing facilities, and is 
not requesting i/o. Thus, the user program will enter 



the i/o library at the point shown by arrow 4 in Figure 
2. Arrows 4-7 show the linkages that will occur. 

Note that each doubleword in hold contains a blank. 
It does not matter whether the blank is to the right or 
to the left, since Fortran data conversion will treat it 
as a zero. (Though if the numbers were whole num- 
bers, it would matter. ) 

To begin with, the user program stores into the 
chcib work area, at chcib2, the address of a parameter 
list .which substitutes for the first four pointers of the 
dcb prefix. 



LA 

L 

ST 



2JPTRS 
3,VCONl 

2,0(0,3) 



Next, since the user program is bypassing chcia, it 
must store the address of the format character string 
into the first word of chcifw and zero out the second 
and third words of chcifw. 

LA 2,FORMAT 

L 3,VCON2 

ST 2,0(0,3) 

SR 4,4 

SR 5,5 

STM 4,5,4(3) 

Then comes the usual sequence of code for calling 

CHCIE. 





LA 


0,200(0,0) 


Indicates that the list item 
is an array, and that the 
array is 200 bytes in 
length. 


CALL < 


CHCIE1, 


(CITEM, 


Causes the conversion and 






ARRAY) 


movement of data to be 
completed. 


HOLD 


DS 


400X 




ARRAY 


DS 


50F 


List Item 


CITEM 


DC 


X'33' 


The first four bits of this 



PTRS 



DC 

DC 
DC 
DC 



A(HOLD) 
A (HOLD) 

A(HOLD-f 

400) 
X'80C00000' 



control byte indicate that 
the elements of the array 
are four bytes long. The 
second four bits indicate 
that the data in the buffer 
is to be converted to real 
form. 

Starting location of raw 
data. 

Current location. ( Same as 
starting location. ) 

End of raw data. 

First byte indicates a 
READ operation. The 
second byte indicates a 
FORMAT statement 
with a list with the 
FORMAT statement 
not encoded. 



VCON1 DC 
VCON2 DC 
FORMAT DC 



V(CHCIB2) 
V( CHCIFW) 
C'(50F8.3)' 
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Appendix C: FORTRAN Data Management 



This appendix describes the assumptions that the For- 
tran i/o library makes in initializing dcbs with infor- 
mation concerning record format (recfm) and data 
set organization (dsorg). These assumptions are de- 
scribed to help reduce a frequent source of error en- 
countered by the user when performing i/o. 

Some introductory material is presented on the dcb 
describing its general use, contents, and sources of 
initialization, before discussing the permissible record 
formats and data set organizations. 



DCB Use 

The Data Control Block (dcb) is created by dcb Man 
agement (chcib) and is used by certain data manage- 
ment routines invoked by macro instruction references 
in i/o Control (chcic). The dcb is required for all i/o 
performed using either bsam or vam. However, the dcb 
is not required for i/o when using the gate macro in- 
structions. 



DCB Content 

The dcb contains information such as the ddname, type 

of data set organization, the type and size of records, 

block size for blocked data sets, number of buffer 

. r-as, exits for synad and eodad, and various control 



;s used by data management. 



DCB Initialization 

The Fortran i/o routines, when processing an input 
data set, take advantage of information in the dcb to 
adapt to the characteristics of the data set and read it 
correctly. Characteristics are based on the parameters 
for a dcb that can be supplied from: 

• The user program — type of i/o used and associated 
c ,ta format. 

• User-supplied ddef commands — some of the infor- 
mation in the dcb can be changed in this manner; 
however, the extent of change is limited. 

• Input data set labels — these override both of the 
above sources of information, within limits set by 
data management. 



Combinations of DSORG and RECFM 

Table 12 gives the permissible combinations of record 
formats and data set organizations that may be speci- 
fied when using the Fortran i/o library. 

Table 12. Combinations of DSORG and RECFM Values 







DSORG VALUES 




RECFM 












vs 


PS 


VSP VI 


VIP 


V 


A 


A 


A 


A 


A 


VB 


N 


A 


N 


N 


N 


VT 


N 


A 


N 


N 


N 


F 


A 


A 


A 


A 


A 


FB 


N 


A 


N 


N 


N 


FS 


N 


A 


N 


N 


N 


FT 


N 


A 


N 


N 


N 


FBS 


N 


A 


N 


N 


N 


FBT 


N 


A 


N 


N 


N 


FBST 


N 


A 


N 


N 


N 


FST 


N 


A 


N 


N 


K 


U 


L 


A 


L 


N 


N 


Codes mean 










A — Acceptable 








L — Limited Acceptable 






N — Not 


acceptable 








DSORG abbreviations mean: 






VS — 


- Virtual sequential ( 


direct-access only ) 




PS - 


- Physical 


sequential 


— BSAM — ( any device 




except te 


rminals ) 






VSP - 


- Virtual sequential partitioned ( like VS ) 




VI - 


- Virtual index sequential ( like VS ) 




VIP - 


- Virtue! index sequential partitioned (like VS ) 


RECFM abb 


reviations 


mean: 






V — 


- Variable- 


length unblocked records 




VB — 


- Variable- 


length blocked records 




VT - 


- Variable-length un 


blocked with track overflow 


F - 


- Fixed-length unblocked records 




FB - 


- Fixed-length blocked records 




FS - 


- Same as 


F, no truncated blocks or unfilled tracks 


FT — 


-Same as 


F, track overflow 




FBS - 


- Same as ! 


7 B, no truncated blocks or unfilled tracks 


FBT — 


- Same as 


FB, track 


overflow 




FBST - 


- Same as 


FBS, tract 


. overflow 




FST -_ 


- Same as 


F, no truncated blocks, track overflow 


U — 


- Undefined record length 





Any of the recfm codes shown can be followed by 
the letter A or M. A indicates that the first character of 
every logical record is an extended axsi fortrax iv 
carriage or punch control code. M indicates that the 
first character of every record is a tss/360 machine con- 
trol byte. In general, the M option cannot be used bv 
Fortran output data sets, since the control codes arc 
unprintable and do not conform to fortrax conven- 
tions. 
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Unformatted FORTRAN Logical Records 

Under any of the organization techniques used, an un- 
formatted write statement may lead to a logical record 
that exceeds the length of the maximum record sup- 
ported by the access method. Furthermore, it is not 
possible to enter the byte size of the entire Fortran 
logical record into the beginning of the i/o physical 
record without the possibility of tying up an indefinite 
amount of virtual storage. Therefore, unformatted 
Fortran logical records may span over data manage- 
ment physical records. In the management of unfor- 
matted Fortran data, the first two bits of every vs 



physical record or the third byte of every ps physical 
record is a control byte defined as follows; 

XW A FORTRAN logical record does not span into or out of 
the data management physical record. 

X'01* This data management physical record is the first of a 
span. 

X'02* This data management physical record is the last of a 
span. 

X*03* This data management physical record is within the 
range of a span. 

No data management physical record will be written 
containing more than one unformatted Fortran logical 
record. 
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Appendix D: DUMP and PDUMP Sample Storage Printouts 



This appendix contains a sample printout for each order: hexadecimal* logical *1, logical *4, integer *2, 

dump format that can be specified for the dump and integer *4, real *4, real *8, complex *8, complex *!$ 

pdump subprogram. The printouts are given in the and literal. 

Table 13. Sample Storage Printouts 



rnNyFR<;TON rnnF n - MFXAr>Fr|MA| 


nomFpp.n cic?c*>cu rscf>C7rR rgmnpn^ r>4D«>Df>D7 


D8D9F7E3 






CONVFRSTON COO* 1 - 1 OGTCAI * 1 


0003E1R0 T F T F T F 


F F 


F F 




CONVERSION CODE 2 - LOGICAL * 4 


0003E1D0 T F T F T F 


CONVERSION CODE 3 - INTEGER * 2 


0003E1RA 1 2 3 


4 


b 6 




CONVERSION CODE 4 - INTEGER * 4 






0003E1F8 12 3 


4 


5 6 


i 

, CONVERSION CODE 5 - REAL * 4 


i 

. 0903E248 1.00000E 00 0.20000E 01 0, 


,30000^ 01 


0.40000E 01 


• 


1 CONVERSION CODE 6 - REAL * 8 

i 







0O03E270 



1.0000CD 00 0.20000D 01 0.30000D 01 0.40000D 01 



IONVERSION CODE 7 - COMPLEX * 8 



O003E2CO 



1.0000CE 00 0.20000E01 Q.20000F Ul 0«30000E 01 



CONVERSION CODE 8 - COMPLEX * 16 



OOO3E310 



1.C0000D 00 0.20000D 01 0.20000H 01 0.30000D 01 



CONVERSION CODE 9 - LITERAL 



C 3E220 ABCDEFGHI JKLM.NOPQRSTUV W 
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Appendix E: Interruption Procedures 



This appendix contains descriptions of the procedures 
followed when the user's program is temporarily inter- 
rupted due to certain types of interrupts. Interrupts 
are hardware-originated breaks in the flow of process- 
ing. Program interrupts result from improper specifica- 
tion or use of instructions and data. The term exception 
is used to refer to these types of interrupts ( see Princi- 
ples of Operation) . Six such exceptions occur* fre- 
quently enough during normal Fortran programming 
to warrant special treatment. 

1. Fixed point overflow exception 

2. Significance exception 

3. Exponent overflow exception 

4. Exponent underflow exception 

5. Floating point-divide exception 

6. Specification exception 

The procedure for handling these exceptions follows. 
The compiler generates code at the beginning of all 
main programs that calls the ciicbdi entry to module 
chcbd. At chcbdi these operations are performed: 

1. Bits are set in the PSW such that the fixed point 
overflow and significance exceptions will be ignored. 

2. Initialization is performed such that control will be 
passed to an entry in module chcbd or chcbe if any 
of the remaining four exceptions occur; 



EXCi ' 


"HON 


ENTRY 


,<-.•!. 1/ 


t.M <How 


CHCBDJ 


>o; "in if 


underflow 


CHCBD4 




. livide 


CHCBD5 


' pCV.ul 


■- .Lion 


CHCBE 1 



At the first three of these entries, flags are set for 
later interrogation by routines called as a result of the 
call overfl ( tests for exponent overflow or under- 
flow, exceptions ) and call dvchk ( tests for floating 
point divide exception ) statements. 

A specification exception occurs when a variable is 
not on a proper word boundary. This condition may 
exist in a Fortran program if an equivalence or com- 
mon statement forces misalignment. The compiler is- 
sues a warning diagnostic, but such a misalignment 
does not prevent the user from executing the program. 
An installation option specifies that one of two courses 
of action is to be taken if a specification interrupt oc- 
curs: (1) terminate the task, or (2) transfer control to 
a program that will perform the desired operation, 
using instructions that will not cause an exception due 
to incorrect boundary alignment. The routine entered 
for either of these eventualities is chcbe, which is 
entered by the chcbei entry. The installation option is 
tested, and one of the two above courses of action 
taken. 

An exponent overflow exception is recognized when 
the result of a floating point addition, subtraction, mul- 
tiplication, or division is either greater than or equal to 
16 63 (approximately 7.2 x 10 75 ). An exponent under- 
flow exception is recognized when the result of a float- 
ing point addition, subtraction, multiplication, or divi- 
sion is less than 16~~ 65 (approximately 5.4 x 10 -79 ). A 
divide exception is recognized when division by zero 
is attempted. 
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Appendix F: Algorithms 



Information about the computations used in the explicitly called mathematical 
subprograms is arranged alphabetically in this appendix, according to subprogram 
module name. The user entry name associated with each subprogram is given in 
parentheses following the module name. 

The information for each subprogram is divided into two parts: a description of 
the algorithm used and a description of the effect of an argument error upon the 
accuracy of the answer ( function value ) . 

The presentation of each algorithm is divided into its major computational 
steps; the formulas necessary for each step are supplied. Some formulas are widely 
known; others are derived from common formulas. In these cases, the process 
leading from the common formula to the computational formula is sketched in 
enough detail that the derivation may be reconstructed. 1 

For the sake of brevity, the needed constants are normally given only symboli- 
cally. (The 1 actual values can be found in the assembly listing of the subprograms. ) 
Some of the formulas are widely known; those that are not so widely known are 
derived from more common formulas. The process leading from the common 
formula to the computational formula is sketched in enough detail so that the 
derivation can be reconstructed by anyone who has an understanding of college 
mathematics and access to the common texts on numerical analysis. 1 Many approxi- 
mations were derived by the so-called "minimax" methods. The approximation 
sought by these methods can be characterized as follows. Given a function f(x), an 
interval I, the form of the approximation (such as the rational form with specified 
degrees), and the type of error to be minimized (such as the relative error), there 
is normally a unique approximation to f(x) whose* maximum error over I is the 
smallest among all possible approximations of the given form. Details of the theory 
and the various methods of deriving such approximation are provided in the 
reference. 1 The accuracy figures cited in the algorithm sections are theoretical, and 
thev do not take round-off errors into account, .vlmor programming techniques 
used to minimize round-off errors are not necessarily described here. 

The accuracy of an answer produced by these algorithms is influenced by two 
factors: the performance of the subprogram and the accuracy of the argument. 
(Performance statistics are given in Table 1.) The effect of an argument error 
upon the accuracy of an answer depends solely upon the mathematical function 
involved and not upon the particular coding used in the subprogram. 

Because argument errors, whether accumulated prior to use of the subprogram 
or introduced by newly converted data, always influence the accuracy of answers, 
a guide to the propagational effect of argument errors is provided. This guide 
(expressed as a simple formula, where possible) is intended to assist users in 
assessing the effect of an argument error. 



1 Any of the common numerical analysis texts may be used as a reference. One such text is F. B. Hildebrand's 
Introduction to Numerical Analysis (McGraw-Hill Book Company, Inc., New York, N.Y., 1956). Background 
information for algorithms that use continued fractions may be found in H. S. Wall's Analytic Theory of 
Continued Fractions (D. VanNostrand C>., Inc., Princc-tim, N.J., 1948). 
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These symbols are used in this appendix to describe the effect of an argument 
error upon the accuracy of the answer: 



Relative error of result given by 
subprogram 



SYMBOL 


EXPLANATION 


gW 


Result given by subprogram 


f(x) 


Correct result 




f(x)-g(x) 


€ 


f(x) 


6 


Relative error of argument 


E 


fW-g(x) 



Absolute error of result given 
by subprogram 



A Absolute error of argument 



The notation used for the continued fractions in this appendix complies with 
the specifications set by the National Bureau of Standards. For more information, 
see Milton Abramowitz and Irene A. Stegun ( editors ) , Handbook of Mathemati- 
cal Functions, Applied Mathematics Series-55 (National Bureau of Standards, 
Washington, D.C., 1965). 

Although it is not specifically stated below for each subroutine, the algorithms 
in this chapter were programmed to conform to the following standards governing 
floating-point overflow/underflow. 

1. Intermediate underflow and overflows are not permitted to occur. This prevents 
the printing of irrelevant messages. 

2. Those arguments for which the answer can overflow are excluded from the 
permitted range of the subroutine. This rule does not apply to cdabs and cabs. 

3. When the magnitude of the answer is less than 16~ 65 , zero is given as the 
answer. If the floating-point underflow exception mask is on at the time, the 
underflow message will be printed. 



Control of Program Exceptions in Mathematical Functions 

The Fortran mathematical functions have been coded with careful control of error 
.!. nations. A rcsul: is provided whenever the answer is within the range repre- 
ss j'iK * iv me iloatrng-p ^n form, .n c. >er to be consistent with Fortran control 

r c? ,; h.c;i[ nver'L'w/undcrfi iw exceptions, tie following types of conditions are 
n .^ OL n z ,w t ; id La; idled sepaiately. 

"vVliv . ; r "»ch r T*ilade of the furcin \. alee is too large to be represented in the 
n .' Jui-ii., thi rendition is c lied a terminal overflow; when the magni- 

fy; . . » 41 > ' co !)c represented, a terminal underflow. On the other hand, if the 

IV.jK-i h;e i reprr: f itablt, but if execution of the chosen algorithm causes an 

; ,-, - r-de ! "Hov t he process, tliir ?oiicIilion is called an intermediate over- 

flow or underflow 

Function subroutines in the Fortran library have been coded to observe the 
following r '-^s for these conditions: 

1. Algorithms which can cause an intermediate overflow have been avoided. 
Therefore an intermediate overflow should not occur during the execution 
of a function subroutine of the library. 

2. Intermediate underflows are detected and not allowed to cause an interrupt. 
In other words, spurious underflow signals are not allowed to be given. Com- 
putation of the function value is successfully carried out. 

3. Terminal overflow conditions are screened out by the subroutine. The argu- 
ment is considered out of range for computation and an error diagnostic is 
given. 
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4. Terminal underflow conditions are handled by forcing a floating-point under- 
flow exception. This provides for the detection of underflow in the same 
manner as for an arithmetic statement. Terminal underflows can occur in the 
following function subroutines: exp, dexp, atan2, datan2, erfc, and derfc. 
For implicit arithmetic subroutines, these rules do not apply. In this case, both 
terminal overflows and terminal underflows will cause respective floating-point 
exceptions. In addition, in case of complex arithmetic (implicit multiply and di- 
vide), premature overflow/underflow is possible when the result of arithmetic is 
very close to an overflow or underflow condition. 
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Explicitly Called Subprograms 
Absolute Value Subprograms 

CABS/CDABS 

1. Write | jc + iy\ = a + ib, 

2. Let v x = max ( jx|, |y| ), and v 2 = min ( |x|, |t/| ). 

3. If characteristics of v^ and v 2 differ by 7 (15 for cdabs) or more, or if v 2 = 0, 
then a = v u b = 0. 

4. Otherwise, 



t)i • Jy 4 + ¥4 (^) 5 and fo = 0. 



If the answer is greater than 16 6S , the floating-point overflow interruption will 
take place (see Appendix C). The algorithms for both complex absolute value 
subprograms are identical Each subprogram uses the appropriate real square root 
subprogram ( sqrt or dsqrt ) . 

Effect of an Argument Error 

1 
€-"2-8. 

Arcsine and Arccosine Subprograms 

ARSIN/ARCOS 

Algorithm 

1. If ^ x ^ V2, then compute arcsin (x) by a continued fraction of the form: 

arcsin (x) = x + x* • F where 

_ d\ d 2 

~(x 2 + d) + (x 2 + c 2 ) ' 

The coefficients of this formula were derived by transforming the minimax 

rational approximation ( in relative error, over the range rg x 2 ^ % ) for 

arcsin (x)/x of the following form; 

arcsin(x) [~ a x + a 2 x 2 "1 

- -a + x~-|^ + feiJc2 + x4 J. 

Minimax was taken under the constraint that a {) = 1 exactly. The relative error 
of this approximation is less than 2~ 2H : \ 
If ^ x ^ %, arccos(x) is computed as: 

arccos ( x ) = — — arcsin ( x ) . 

2. If V2 < x^ 1, then compute arccos(x) essentially as: 

arccos(x) = 2 • arcsin I .* / — - — I. 
This case is now reduced to the first case because within these limits, 



1 - x 

This computation uses the real square root subprogram (sqrt) 
If V2 < x ^ 1, arcsin(x) is computed as: 

arcsin ( x ) = -— — arccos ( x ) . 

Implementation of the above algorithms ( steps 1 and 2 ) were carried out with 
care to minimize the round-off errors. 
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3. If — 1 fi .t < 0, then arcsin ( x ) = — arcsin jxj 
andarecos(x) = ir — arccos \x\. 
This reduces these cases to one of the two positive eases. 



Effect of an Argument Error 

A 



i. For small values of x, E ~ A. Toward the limits ( ± 1) of the 



V 1 - x- 

range, a small A causes a substantial error in the answer. For the arcsine, e ~ o 
if the value of x is small, . 



DARSIN/DARCOS 

Algorithm 

1. If :g x g %, then compute arcsin(x) by a continued fraction of the form: 

arcsin (x) ^x + x 3, F where 

d\ d 2 d$ d± 



F = c x + 



(x 2 + c 2 ) + (x 2 + c H ) 4- (x- + c 4 ) + (x 2 + c 5 )" 



The relative error of this approximation is less than 2~ 57 - 2 . 

The coefficients of this formula were derived by transforming the minimax 
rational approximation ( in relative error, over the range rg x 2 ^ % ) for 
arcsin ( x ) / x of the following form : 

arcsin(x) S 0,1 "*" flaX " "*" a * xA "*" a ^ ~*~ fli > x *~] 

- - o + x^ ^ + fe ^ 3 + h ^ + fe ^ + x8 J . 

Minimax was taken under the constraint that a i} — 1 exactly. 
If :g x :g %, arccos(x) is computed as: 

arccos(x) = — arcsin(x). 

2. If % < x rg 1, then compute arccos (x) essentially as: 
arccos(x) = 2 • arcsin 1 ^\ — - — I. 
This case is now reduced to the first case because within these limits, 



<>aj 



1 -x 



2 = 



^y 2 . 



This computation uses the real square root subprogram (dsqrt). 
If V2 < x ^ 1, arcsin(x) is computed as: 

arcsin(x) = — arccos(x). 

Implementation of the above algorithms (steps 1 and 2) were carried out with 
care to minimize the round-off errors. 
3. If — 1 :g x < 0, then arcsin(x) = — arcsin |x|, and arccos(x) = tt — arccos |xl. 
This reduces these cases to one of the two positive cases. 

Effect of an Argument Error 
A 

^ ~ /] — .-j * F° r small values of x, E ~ A. Toward the limits ( ± 1 ) of the 

range a small A causes a substantial error in the answer. For the arcsine, c ~ 8 if 
the value of x is small. 
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Arctangent Subprograms 
ATAN/ATAN2 

Algorithm 

1. For arctan (x u x 2 ): 

If Xx < 0, use the identity aretan (x u x 2 ) = —arctan ( — Xi, x 2 ). 
Hence we may assume that x { ^ 0. Then: 



If either x 2 = or 

T Xi 
If x 2 < and 



x-> 



> 2 24 , the answer = — 



< 2~~ 24 , the answer = ?r. 

For the general case, if x 2 > 0, the answer = arctan I FH 1 , and 
if x 2 < 0, the answer = ?r — arctan I h—j I . 

The remainder of the computation is identical for either one or two arguments. 

2. Reduce the computation of arctan (x) to the case rg x :g 1, by using 

arctan ( — x) = — arctan (x), or 

arctan I t~t I =-^ — arctan !x|. 
\[x|/ 2 

3. If necessary, reduce the computation further to the case |x| 5^ tan 15° by using 

/y~3-x- 1\ 

arctan (x) = 30° + arctan I ; -=- 1. 

_ \ x + V3 / 

\/3-x - 1 



The value of , 

x + \f 3 

tan 15° < x 5i 1. The va 



tan 15° if the value of x is within the range, 
ue of ( V 3 • x — 1 ) is computed as 



( V 3 — l) x — l + xto avoid the loss of significant digits. 

4. For |x[ ^ tan 15°, use the approximation formula: 

arctan ( x ) 55913709 
— = 0.60310579 - 0.05160454x 2 4- v° D * lotx ™ 



x 2 + 1.4087812 * 



This formula has a relative error less than 2~~ 27 - 1 and can be obtained by 
transforming the continued fraction 

arctan (x) x -» 5 



x 3 + 

w 



(*— )- 

— -==x~ 2 H — — — I 10~ 4 , but the true 



4-5 

, , . 7-7-9 

value or ic is 



(rSi + -) 



+ x- 2 I + 



The original continued fraction can be obtained by transforming the Taylor 
series into continued fraction form. 



42 



Effect of an Argument Error 

E ~ -r-r — 7. For small values of x, c ~ S; as the value of x increases, the effect 

1 + X" 

of 8 upon c diminishes. 

DATAN/DATAN2 

Algorithm 

L For arctan (xi, x 2 ): 

If x x < 0, use the identity arctan (xi, x 2 ) = — arctan ( — x u x 2 ). 
Hence we may assume that x a ^ 0. Then: 



If either x 2 = or 



> 2 50 , the answer = — . 



If x 2 < and 



*1 



*2 



< 2~~ 56 , the answer = w: 



For the general case, if x 2 > 0, the answer = arctan I p-j I , and 
if x 2 < 0, the answer = * — arctan 1 p-l j. 

The remainder of the computation is identical for either one or two arguments. 

2. Reduce the computation of arctan (x) to the case ^ x ^ 1 by using 

arctan ( — x) = — arctan(x) and 

arctan , , =-r — arctan xL 

|x| 2 

3. If necessary, reduce the computation further to the case |x| ^ tan 15° by using 

, , /Y3-X- 1\ 
arctan(x) = 30° + arctan i ; j=- I. 

V3 « x - 1 
The value of ; j=~ < tan 15°, if the value of x is within the range tan 

x + Y3 ~~~ 

15° < x ^ 1. The value of ( V3 • x - 1) is computed as ( Y~3 - 1) x - 1 + x 

to avoid the loss of significant digits. 

The relative error of this approximation is less than 2 -60 7 . 

The coefficients of this formula were derived by transforming a minimax 
rational approximation ( in relative error, over the range ^g x 2 ^ 0.071797 ) 
for arctan(x)/x of the following form: 

arctan(x) [~c +C1X 2 + c 2 x* + c 3 x°l 

- = oo + x- L do + djx 2 + d 2 x^ +x«J 

Minimax was taken under the constraint that a {) = 1 exactly. 

4. For |x| ^ tan 15°, use a continued fraction of the form: 

arctan(x) ^ F Oi <h Q 3 1 
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Effect of an Argument Error 

A 



^ - .,. For small values of x, e ~ 8, and as the value of x increases, the effect 

of c upon 8 diminishes. 

Error Functions Subprograms 

ERF/ERFC 

Algorithm 

1. If 5^ x g 1, then compute the error function by the following approximation: 

erf(x) e= x(a -f a x x 2 + a 2 * 4 + . . . + fls* 10 ). 

The coefficients were obtained by the minimax approximation (in relative 
error) of erf(x)/x as a function of x 2 over the range ^ x 2 ^ 1. The relative 
error of this approximation is less than 2~ 24 - 6 . The value of the complemented 
error function is computed as erf c(x) = 1 — erf(x). 

2. If 1 < x < 2.040452, then compute the complemented error function by the 
following approximation: 

erfc (x) = foo + b x z + b 2 z 2 + . . . + & 9 z 9 

where z = x — T and T () ~ 1.709472. The coefficients were obtained by the 
minimax approximation (in absolute error) of the function f(z) = erfc (z + 
T ) over the range - 0.709472 ^ z ^ 0.33098. The absolute error of this 
approximation is less than 2~ 31 - 5 . The limits of this range and the value of the 
origin T were chosen to minimize the hexadecimal round-off errors. The value 

1 

of the complemented error function within this range is between -r^ and 0.1573. 

The value of the error function is computed as erf (x) = 1 — erfc(x). 

3. If 2.040452 ^ x < 13.306, then compute the complemented error function by 
the following approximation: 

erfc(x) = e _r • Fix where z = x 2 and 
C\ ■: c 2 z + c- d z 2 



F = co + 



cUz + cLz 2 + 



The coefficients for F were obtained by transforming a minimax rational 
approximation (in absolute errors, over the range 13.306 -2 ^ w ^ 2.040452 -2 ) 

of the function f(w ) = erfc( x ) • x • e /£ , to = x~ 2 , of the following form: 



a {) + a x w + a>w 2 + a 3 ti; 3 

f(io) =■ 



b + biiv 4- iv 2 



The absolute error of this approximation is less than 2~ 261 . This computation 
uses the real exponential subprogram ( exp ) . 
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If 2.040452 ^ x < 3.919206, then the error function is computed as 

erf(x) = 1 — erfc(x). 

If 3.919206 g x, then the error function is ~ 1. 

4. If 13.306 fg x, then the error function is a= 1, and the complemented error func- 
tion is a* ( underflow ) . 

5. If x < 0, then reduce to a case involving a positive argument by the use of the 
following formulas: 

erf( — x) = — erf(x), anderfc ( — x) = 2 — erfc(x). 

Effect of an Argument Error 

E ~ e ~ x * • A . For the error function, as the magnitude of the argument exceeds 1, 
the effect of an argument error upon the final accuracy diminishes rapidly. For 
small values of x, c ~ 8. For the complemented error function, if the value of x is 

-. Therefore, c ~ 2 x 2 • 8. If the value of x is negative 



greater than 1, 


erfc(x) 


*-*•» 


2x 


or less than 1, then « ~ 


e - 


-* 2 • A. 


DERF/DERFC 









Algorithm 

1. If :g x < 1, then compute the error function by the following approximation: 

erf(x) as x(a + a x x 2 + a 2 x* -f . . . + fl n x 22 ). 

The coefficients were obtained by the minimax approximation (in relative 
error) of erf (x)/x as a function of x 2 over the range ^ x 2 ^ 1. The relative 
error of this approximation is less than 2~ 56 9 . The value of the complemented 
error function is computed as erfc(x) = 1 — erf(x). 

2. If 1 ^ x < 2.040452, then compute the complemented error function by the 
following approximation: 

erfc(x) = fo + b x z -f b 2 z 2 + . . . fc 18 z 18 

where z = x - T and T ~ 1.709472. The coefficients were obtained by the 
minimax approximation (in absolute error) of the function f(z) = erfc(z + T ) 
over the range -0.709472 ^ z ^ 0.33098. The absolute error of this approxi- 
mation is less than 2~ 60 - 3 . The limits of this range and the value of the origin 
T were chosen to minimize the hexadecimal round-off errors. The value of the 

complemented error function within this range is between -^r- and 0.1573. The 

Zoo 

value of the error function is computed as erf (x) = 1 — erfc(x). 

3. If 2.040452 ^ x < 13.306, then compute the complemented error function by 
the following approximation: 

erfc(x) a^ e~ z • F/x where z = x 2 and 

p = , d x d 2 d$ d 7 

°° (z + c,)+ (z + c 2 ) + "* (z + c 6 ) + (z + c 7 )' 

The coefficients for F were derived by transforming a minimax rational approxi- 
mation (in absolute errors, over the range 13.306" 2 ^ w ^ 2.040452" 2 ) of the 
function f(w) = erfc(x) • x • e **, w = x" 2 , of the following form: 

a -f a x w -f a 2 w 2 + . . . + a 7 w 7 



f(w) a. 



b# + b x w + b 2 w 2 + . . . + betv 6 +W 1 ' 



The absolute error of this approximation is less than 2~ 57 9 . This computation 
uses the real exponential subprogram (dexp). If 2.040452 ^ x < 6.092368, then 
the error function is computed as erf (x) = 1 — erfc(x). 
If 6.092368 g x, then the error function is ~ 1. 
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4. If 13.306 g x, then the error function is =s 1, and the complemented error 
function = (underflow). 

5. If x < 0, then reduce to a case involving a positive argument by the use of the 
following formulas; 

erf( — x) = — erf(x), and erfc ( — x) = 2 — erfc(x). 

Effect of an Argument £ rror 

E ~ e~ r ~ • A. For the error function, as the magnitude of the argument exceeds 
1, the effect of an argument error upon the final accuracy diminishes rapidly. For 
small values of x, e ~ 8. For the complemented error function, if the value of 

e -x2 

x is greater than 1, erfc(x) ~ — - — . Therefore, c ~ 2x 2 • 8. If the value of x 
is negative or less than 1, then € ~ e~ x2 • A. 



Exponential Subprograms 
EXP 

Algorithm 

1. If x < — 180.218, then is given as the answer via floating-point underflow. 

2. Otherwise, divide x by log ( >2 and write 

x 

V = l n = 4a — b — d 

log 2 

where a and & are integers, 0^6g3 and ^ d < 1. 

3. Compute 2~ (i by the following fractional approximation; 

2d 



1 - 



617.97227 
0.034657359 d 2 + d + 9.9545948 - - 



d? 4- 87.417497 
This formula can be obtained by transforming the Gaussian continued fraction 

zzzzzzzz 



1+ 2- 3+ 2- 5+ 2- 7+ 2* 

The maximum relative error of this approximation is 2~~ 29 . 

4. Multiply 2~ rf by 2-'\ 

5. Finally, add the hexadecimal exponent a to the characteristic of the answer. 

Effect of an Argument Error 

c — J If the magnitude of x is large, even the round-off error of the argument 
causes a substantial relative error in the answer because A = $ • x. 



DEXP 

Algorithm 

1. If x < — 180.2187, then is given as the answer via floating-point underflow. 

2. Divide x by Iog e 2 and write 

x = I 4 a - b - -~\ • log c 2 - r 
where a, 6, and c are integers, g Z? ^ 3, g c g 15, and the remainder r is 
within the range g r < — • log e 2. This reduction is carried out in an extra 
precision to ensure accuracy. Then e**' = 16 a • 2~ b • 2~ c/16 • e~ r . 
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3. Compute e~ r by using a minimax polynomial approximation of degree 6 over 

the range fg r < -75- • loge2. In obtaining coefficients of this approximation, 
16 

the minimax of relative errors was taken under the constraint that the constant 

term Oo shall be exactly 1. The relative error is less than 2 -56 * 87 . 

4. Multiply e~ r by 2~ c/ie . The 16 values of 2~ e/16 for ^ c 5g 15 are included in 
the subprogram. Then halve the result b times. 

5. Finally, add the hexdecimal exponent of xr to the characteristic of the answer. 

Effect of an Argument Error 

E ~ A. If the magnitude of x is large, even the round-off error of the argument 
causes a substantial relative error in the answer because A = 8 • x. 

CEXP/CDEXP 
Algorithm 

The value of e x + iy is computed as e x • cos(y) + i • e x • sin(y). The algorithms for 
both complex exponential subprograms are identical. Each subprogram uses the 
appropriate real exponential subprogram ( exp or dexp ) and the appropriate real 
sine/ cosine subprogram (cos/ sin or dcos/dsin). 

Effect of an Argument Error 

The effect of an argument error depends upon the accuracy of the individual 
parts of the argument. If e*+^ = R * e iH 7 then H = y and c(R) ~ A (x). 



Gamma and Log Gamma Subprograms 

GAMMA/ ALGAMA 
Algorithm 

1. If < x ^ 2~ 252 , then compute log-gamma as log e r(x) ss — log e (x). 
This computation uses the real logarithm subprogram ( alog). 

2. If 2~ 252 < x < 8, then compute log-gamma by taking the natural logarithm of 
the value obtained for gamma. The computation of gamma depends upon the 
range into which the argument falls. 

T(x + 1) 

3. If 2~ 252 < x < 1, then use T(x) = to reduce to the next case. 

x 

4. If 1 :g x g 2, then compute gamma by the minimax rational approximation ( in 
absolute error) of the following form; 

. . z [a + atf + a 2 z 2 + azz z ] 

T(x) » c + — , — 

Z?o + fciz + O22 2 + z 3 

where z = x — 1.5. The absolute error of this approximation is less than 2 -25 - 9 . 

5. If 2 < x < 8, then use T(x) = (x— l)r(x-l)to reduce step by step to the 
preceding case. 

6. If 8 ^ x 9 then compute log-gamma by the use of Stirling's formula: 

log e r(x) -x(log e (x) - 1) - %log e (x) + %log e (2 ? r) + G(x). 
The modifier term G(x) is computed as 
G(x)^d x- 1 + d!X- 2 . 

These coefficients were obtained by a form of minimax approximation minimiz- 
ing the ratio of the absolute error to the value of x. The absolute error is less 
than x • 2~ 26 - 2 . Remembering the fact that x < log e r(x) in this range, the 
con! ribution of this error to the relative error of the value for log-gamma is less 
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than 2~~ 2G - 2 . This computation uses the real logarithm subprogram ( alog). 

For gamma, compute r(x) = e v > where y is the value obtained for log-gamma. 

This computation uses the real exponential subprogram (exp). 

Effect of an Argument Error 

c — ' \p (x) • A for gamma, and E ~ xj/ (x) • A for log-gamma, where \p is the 
digamma function. 

If — < x < 3, then -2 < $ (x) < 1. Therefore, E ~ A f or log-gamma. How- 

ever, because x = 1 and x = 2 are zeros of the log-gamma function, even a small 
8 can cause a substantial e in this range. 

If the value of x is large, then \p (x) ~ log e (x). Therefore, for gamma, 
e ~ 8 x • log e (x). In this case, even the round-off error of the argument con- 
tributes greatly to the relative error of the answer. For log-gamma with large 
values of x, c ^ B. 



DGAMMA/DLGAMA 

Algorithm 

1. If < x rg 2~ 252 , then compute log-gamma as log e r(x) ^ — log e (x). 
This computation uses the real logarithm subprogram ( dlog ) . 

2. If 2~ 252 < x < 8, then compute log- gamma by taking the natural logarithm 
of the value obtained for gamma. The computation of gamma depends upon the 
range into which the argument falls. 

r(x + 1) 
3 jf 2-252 < x < 1, then use r(x) = to reduce to the next case. 

4. If 1 5^ x rg 2, then compute gamma by the minimax rational approximation (in 
absolute error) of the following form: 

z [gp + a<z + . . . + a Q z G ] 
T(x) = Co + ^ + fciZ + _ + bQz6 + %1 

where z = x — 1.5. The absolute error of this approximation is less than 2 -59 3 . 

5. If 2 < x < 8, then use r(x) = (x — 1) r(x — 1) to reduce to the preceding 
case. 

6. If 8 ^ x, then compute log-gamma by the use of Stirling's formula: 

log e r(x) -X(l0g t ,(x) - 1) - y 2 l0g e (x) + %l0ge(&r) +G(x). 

The modifier term G(x) is computed as 

G(x) ^ d x-~i + dix-* + rf 2 x- 5 + d 3 x~« + d 4 x~ { \ 

These coefficients were obtained by. a form of minimax approximation minimiz- 
ing the ratio of the absolute error to the value of x. The absolute error is less 
than x • 2 -56 - 1 . Remembering the fact that x < log e r(x) in this range, the 
contribution of this error to the relative error of the value for log-gamma is less 
than 2 -56 - 1 . This computation uses the real logarithm subprogram (dlog). For 
gamma, compute r(x) = e ?/ , where y is the value obtained for log-gamma. 
This computation uses the real exponential subprogram (dexp). 

Effect of an Argument Error 

c ~ </f(x) • A for gamma, and E ~ i//(x) • A for log-gamma, where \j/ is the 
digamma function. 

1 

If -y < x < 3, then — 2 < \j/(x) < 1. Therefore, E ~ A f or log-gamma. How- 
ever, because x = 1 and x = 2 are zeros of the log-gamma function, even a small 
8 can cause a substantial e in this range. 

48 



If the value of x is large, then <//(%) ~ log e (x). Therefore, for gamma, 
c ~ 8 • x • loge(x). In this case, even the round-off error of the argument con- 
tributes greatly to the relative error of the answer. For log-gamma with large 
values of x, c ^ 8. 



Hyperbolic Sine and Cosine Subprograms 

S1NH/COSH 

Algorithm 

1. If |z| < 1.0, then compute sinh(x) as: 

sinh(x) sag x + CiX 3 + c 2 x 5 + c 3 x 7 . 

The coefficient c { were obtained by the minimax approximation (in relative 

sinh(x) 
error) of as the function of x 2 . The maximum relative error of this 

1 x 

approximation is 2~ 25 - 8 . 

2. If x ^ 1.0, then sinh(x) is computed as: 

sinh(x) = (1 + 8) [e x+log e v - t> 2 /e*+ l0B . v ]. 

Here, 1 + 8 = — , so that this expression is theoretically equivalent to 

ZtV) 

[e x — e~ x ]/2. The value of v (and consequently those of log e t> and 8) was so 
chosen as to satisfy the following conditions: 

a ) v is slightly less than %, so that 8 > and small. 

fe ) log e D is an exact multiple of 2~ 16 . 
The condition b) insures that the addition x + log fc t; is carried out exactly. This 
maneuver was designed to reduce the round-off errors and also to enlarge the 
limits of acceptable arguments. This computation uses the real exponential sub- 
program (exp). 

3. If x ^ — 1.0, use sinh( x) = — sinh( |x| ) to reduce to case 2 above. 

4. If cosh (x) is desired, then for all valid values of arguments use the identity: 
cosh(x) = (1 + 8) [e^ 10 *.* + v 2 /e*+ log e v ]. Here the notation and the consid- 
eration are identical to case 2 above. This computation uses the real exponential 
subprogram (exp). 

Effect of an Argument Error 

For the hyperbolic sine, E ~ A • cosh ( x ) and c ~ A • coth ( x ) . 
For the hyperbolic cosine, E ~ A • sinh(x) and c ~ 8 • tanh(x). 

Specifically, for the cosine, c ~ A over the entire range; for the sine, c ~ 8 f or 
small values of x. 



DSINH/DCOSH 

Algorithm 

1. If |x| < 0.881374, then compute sinh(x) as: 

sinh (x) ==£CqX + dx 3 + c 2 x 5 + . . . + c 6 x 13 . 

The coefficients c { were obtained by the minimax approximation (in relative 
sinh(x) 

error ) of as the function of x 2 . Minimax was taken under the constraint 

x 

that c = 1 exactly. The maximum relative error of this approximation is 2~ 55 - 7 . 

2. If x ^ 0.881374, then sinh(x) is computed as: 

sinh(x) = (1 + 8) [e x+]os * v - v 2 /e s + los * v ]. 
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Here, 1 + 8 = — , so that this expression is theoretically equivalent to 

[e r — e-^/2. The value of v (and consequently those of log e t; and 8) was so 
chosen as to satisfy the following conditions: 

a ) vis slightly less than %, so that S > and small. 

b ) log e t> is an exact multiple of 2 _ir \ 

The condition b ) insures that the addition x + log e u is carried out exactly. This 
maneuver was designed to reduce the round-off errors and also to enlarge the 
limits of acceptable arguments. This computation uses the real exponential sub- 
program (dexp). 

3. If x fg —0.881374, then use sinh(x) = — sinh( |x| ) to reduce to case 2 above. 

4. If cosh(x )is desired, then, for all valid arguments use the identity: 

eosh(x) = (1 + 8) [e- r + los « r + t; 2 /e a + Iog e u ]. Here the notation and the consid- 
eration are identical to case 2 above. This computation uses the real exponential 
subprogram (dexp). 

Effect of an Argument Error 

For the hyperbolic sine, E ~ A • cosh(x) and e ~ A • coth(x). 
For the hyperbolic cosine, E ~ A • sinh(x) and c ~ A • tanh(x). 

Specifically, for the cosine, e — A over the entire range; for the sine, * ~ 8 for 
the small values of x. 



Hyperbolic Tangent Subprograms 

TANH 

Algorithm 

1. If |x|^ 2- ,L> , then tanh(x) = x. 

2 jf 2~~ 12 < |x| ^ 0.7, use the following fractional approximation: 
tanh(x) f 0.8145651 



f 0.8145651 "I 

l-40.0037828 + ^ + 2-471749 j. 



The coefficients of this approximation were obtained by taking the minimax 
of relative error, over the range x 2 < 0.49, of approximations of this form under 
il? constraint 'hat the first term shall be exactly 1.0. The maximum relative 

error of this approximation is 2~ 26 * 4 . 

2 

.1 If" 0.7 <x< 9.011, then use the identity tanh(x) = 1 - - t — , Q t , . 

y x ; (e x ) 2 + 1 

The :o./.put;;tion for this case uses the real exponential subprogram (exp). 
•i. Vi x _j ?/.()ll, fhantanh(x) ^ 1. 
5. Ifx< —0.7, then use the identity tanh(x) = — tanh( — x). 

Effect of an Argument Error 

2a 
E ~ (1 — tanh 2 x) A, and € ~ , , /rt x . For small values of x, c — ' 8, and as the 

smh (2x) 

value of x increases, the effect of 8 upon e diminishes. 

DTANH 

Algorithm 

1. If |x| ^ 2~- 2K , then tanh(x) = x. 

2 If 2" 2S < jx| < 0.54931, use the following fractional approximation: 
tanh(x) d x x 2 d 2 d s 

= Co H - ; ; - ; ; r ; , 

x x z + d + x- + Co + x 2 + c 3 
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This approximation was obtained by rewriting a minimax approximation of the 
following form : 

tanh(x) a {) + a x x 2 + avx 4 

£* C + X 2 * 



x " " & + ki* 2 + b 2 x* + x 6 ' 

Here the minimax of relative error, over the range x 2 5^ 0.30174, was taken 
under the constraint that c shall be exactly 1.0. The maximum relative error of 
the above is 2~ m . 

2 

3. If 0.54931 <i x < 20,101, then use the identity tanh(x) = 1 - 2x . 

This computation uses the double precision exponential subprogram ( dexp ) . 

4. If x ^ 20.101, then tanh(x) =- 1. 

5. If x^ — 0.54931, then use the identity tanh(x) = - tanh( - x). 

Effect of an Argument f rror 

2 A 
E — (1 - tanh 2 x) A, and c ~ — , ,„ . For small values of x, c ~ 8. As the 

sinh ( 2x ) 

value of x increases, the effect of 8 upon € diminishes. 



Logarithmic Subprograms (Common and Natural) 

ALOG/ALOG10 
Algorithm 

1. Write x = l^ • 2~ r/ • m where p is the exponent, </ is an integer, ^ q ^ 3, 
and m is within the range, % 5^ m < 1. 

2. Define two constants, a and b (where a = base point and 2~ h = a), as follows: 

If % <* m < — = then a = % and b = 1. 

V2 

1 

If ~~= ^ m < 1, then a = 1 and & = 0. 
y A 

rt „ T m — a m 1 + z 

3. Write z = — — -. Then, m = a- - and \z\ < 0.1716. 

m + a 1 — z 

1 + z / 1 + z\ 

4. Now,x = 2 4 ^-^-^ j— — ,andlog e (x) = (4p - qr - fo) lofo 2 + ^g J t-^- 1 

5. To obtain log J Z\ e rs t compute u; = 2z = „ ™ , ** g (which is rep- 

\1 - z/' F 0.5m + 0.5a 

resented in our system with slightly more significant digits than z itself), and 
apply an approximation of the following form: 

1o& - (t^) ^ w [ c » + ;r=y • 

These coefficients were obtained by the minimax rational approximation of 

1 . /I + z\ 
2~log e l y— — I over the range z" 2 c (0, 0.02944) under the constraint that c 

shall be exactly 1.0. The maximum relative error of this approximation is less 
than2~ 25r « 

6. If the common logarithm is desired, then log H ,x = log 10 e • log«Jt. 

Effect of an Argument Error 

E ~ 8. Specifically, if 8 is the round-off error of the argument, e.g., $ ~ 6 • 10~ 8 , 
then E ~ 6 • 10~ 8 . Therefore, if the argument is close to 1, the relative error can 
be very large because the value of the function is very small. 
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DLOG/DLOG10 

Algorithm 

1. Write x = 16 P • 2~ r/ 8 m where p is the exponent, q is an integer, ^ q ^ 3, 
and m is within the range V2 ^ w < 1. 

2. Define two constants, a and fo (where a = base point and 2~ b = a), as follows: 

1 

If y 2 < m < "T^ 1 , then a = % and Zr = 1. 

— V2 

If — = < m < 1, then (7=1 and & = 0. 

V2 — 

m — a 1 + z 

3. Write z = - — . Then, m = a* - and \z\ < 0.1716. 

m -f a 1 — z 

1 + z /l + z\ 

4. Now, x = 2 4p ~ ff - ft • _ , and Io&jc = (4p - q - b) log e 2 + log e I I. 

(1 + z\ m — a 

__ I , first compute tv = 2z = — — fl ^ ( which is repre- 
sented in our system with slightly more significant digits than z itself), and 
apply an approximation of the following form : 

log t . I , _ 1 = w c + du;- / w 1 + c 2 + : — \ . 

V Z/ I io^+c+ + -ri—f 

— \ tv- + C G / -J 

These coefficients were obtained by the minimax rational approximation of 

1 /l + z\ 

— log, I _ ^ I over the range z~ e ( 0, 0.02944 ) under the constraint that c„ 

^z \ I z / 

shall be exactly 1.0. The maximum relative error of this approximation is less 
than2-«°- 5r \ 
6. If the common logarithm is desired, then log H »x = log 10 e • log e x. 

Effect of an Argument Error 

E ~ 8. Therefore, if the value of the argument is close to 1, the relative error can 
be very large because the value of the function is very small. 



CLOG/CDLOG 

Algorithm 

1. Write log,, (x 4- iy) = a + ih. 

2. Then, a = log,. |x -f iy\ and b = the principal value of arctan (y, x). 

3. log,. jx + iy\ is computed as follows: 

Let d, = max ( ! x|, !{/j ), and v 2 — min ( |x|, \y\ ). 

1 

Let t be the exponent of v u i.e., V] = m • 16 f , -r^ ^ jn < 1- 



Finally, let t 



\tiit^0 \ 



it it t ^ u 
\f — 1 if * 



1 - w - nff> of. 
and s = 16^. 



Then, log, |x + iy\ = At x • log,(2) + % log, [(^) 2 + (-7-)" ]• 

Computation of vjs and v^/s are carried out by manipulation of the charac- 
teristics of Vi and v 2 . In particular, if v-J s < 1, it is taken to be 0. The algor- 
ithms for both complex logarithm subprograms are identical. Each subprogram 
uses the appropriate real natural logarithm subprogram ( aloc or dloc) and the 
appropriate arctangent subprogram (atan2 or datan2). 
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Effect of on Argument Error 

The effect of an argument error depends upon the accuracy of the individual 
parts of the argument. If x + iy = r • e** and loge (x + iy) = a + ib, then h = h 
andE(a) = B(r). 



Sine and Cosine Subprograms 
SIN/COS 

Algorithm 

4 

1. Define z = — * |x| and separate z into its integer part (q ) and its fraction part 

(r). Then z = q + r, and |x| = I -j- • q 1 + I -j- • r ■• 

2. If the cosine is desired, add 2 to q. If the sine is desired and if x is negative, 
add 4 to q. This adjustment of q reduces the general case to the computation 
of sin ( x ) for x ^ because 

cos ( ±x) = sin I ~-+ x J, and 
sin ( — x) = sin(ir + x). 

3. Let qo = q mod 8. 

Then, for q Q = 0, sin (x) = sin I — • r I , 

q Q = 1, sin (x) = cos ("^(l " f ) h 

qr = 2, sin (x) = cos (-J-**J> 

q = 3, sin (x) = sin I * (1 - r)Y 

q = 4, sin (x) = - sin f — • r j, 

g = 5, sin (x) = - cos ("^(l - r) h 

q = 6, sin (x) = — cos I — • r L 

q = 7, sin (x) = - sin ("J-(1 - r)J. 
These formulas reduce each case to the computation of either sin 1 — • r x \ 

or cos 1 — • fi I where r x is either r or ( 1 — r ) and is within the range, 
^ n ^ 1. 

4. If sin I -j- • r x lis needed, it is computed by a polynomial of the following 
form: 

sin ("T * n / fl ^ a ° + fllfl2 + a2fl4 + flafl ^* 

The coefficients were obtained by the interpolation at the roots of the Chebyshev 
polynomial of degree 4. The relative error is less than 2 -28 - 1 for the range. 
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5. If cos 1 — • n lis needed, it is computed by a polynomial of the following 
form: 

cos ( -J- • n )= 1 + fo^! 2 + ZW +fe 3 ri G . 

Coefficients were obtained by a variation of the minimax approximation which 
provides a partial rounding for the short precision computation. The absolute 
error of this approximation is less than 2~ 24 - 57 . 

Effect of an Argument Error 

E ~ A. As the value of x increases, A increases. Because the function value dimin- 
ishes periodically, no consistent relative error control can be maintained outside 

the principal range, — — ^ x ^ + — . 
2t ii 



DSIN/DCOS 

Algorithm 

IT 

1. Divide \x\ by - " and separate the quotient (z) into its integer part (q) and 

. 4 
its fraction part (r). Then, z = \x\ • — = q + r , where q is an integer and r 

is within the range, fg r < 1. 

2. If the cosine is desired, add 2 to q. If the sine is desired and if x is negative, 
add 4 to q. This adjustment of q reduces the general case to the computation of 

sin ( x ) for x ^ 0, because 

cos ( ± x) = sin I \x\ + — I, and 

sin ( — x) = sin ( \x\ + tt). 

3. Let f/ == q mod 8. 

Then, for q = 0, sin (x) = sin I — • r I, 

r/o = l,sin (x) = cos I -^- ( 1 - r)j, 
c/o = 2, sin (x) = cos I — • r 1, 
fld = 3, sin (x) = sin I — (1 — r) J, 
q n — 4, sin (x) = — sin I — • r I, 
q i} = 5, sin (x) = - cos ( -^- ( 1 - r)j, 
q {) — 6, sin ( x ) = — cos ( — • r \ , 
q i} = 7, sin (x) = - sin (^~(1 ~ r)\. 
These formulas reduce each case to the computation of either sin 



or cos 

^ n ^ 1 
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I -T- • fi 1 ; where r x is either r or ( 1 — r ), and is within the range, 



(t-0 



4. Finally, either sin I — • r x lor cos I ~T~ * f 1 I is computed, using the polynomial 

interpolations of degree 6 in r x 2 for the sine, and of degree 7 in r x 2 for the cosine. 
In either case, the interpolation points were the roots of the Chebyshev poly- 
nomial of one higher degree. The maximum relative error of the sine polynomial 
is 2 -58 and that of the cosine polynomial is 2 -84 3 . 

Effect of on Argument Error 

E ~ A. As the value of the argument increases, A increases. Because the function 
value diminishes periodically, no consistent relative error control can be main- 
tained outside of the principal range, — — ^ x fg +-<>-• 

CSIN/CCOS 
Algorithm 

1. If the sine is desired, then 

sin(x + it/) = sin(x) • cosh(y) + i 9 cos(x) • sinh(^). 
If the cosine is desired, then 

cos(x + iy) = cos (x) • cosh(y) — i* sin(x) • sinh(y). 

2. The value of sinh(x) is computed within the subprogram as follows. 
Assume x ^ for this, since sinh( — x) = — sinh(x). 

3. If x^ 0.346574, then use sinh (x) = Vz (e* -Y 

4. If fg x < 0.346574, then compute sinh( x) by use of a polynomial; 

sinh(x) 

a* a + a x x 2 + a 2 x*. 

The coefficients were obtained by the minimax approximation (in relative 
error) of sinh(x)/x over the range ^ x 2 ^ 0.12011 under the constraint that 
a shall be exactly 1.0. The relative error of this approximation is less than 

2-26.18 



e l*|* 



5. The value of cosh(x) is computed as cosh (x) = sinh|x| + 

This computation uses the real expoential subprogram (exp) and the real 
sine/ cosine subprogram ( sin /cos ) . 

Effect of on Argument Error 

To understand the effect of an argument error upon the accuracy of the answer, the 
programmer must understand the effect of an argument in the sin/cos, exp, and 
sinh/ cosh subprograms. 



CDSIN/CDCOS 

Algorithm 

1. If the sine is desired, then 

sin (x + iy) = sin(x) • cosh(t/) + i • cos(x) • smh(y). 
If the cosine is desired, then 

cos(x ■+ iy) =x cos(x) »cosh(i/) - %• sin(x) •sinh(y). 

2. The value of sinh(x) is computed within the subprogram as follows. 
Assume x ^ for this, since sinh( — x) = - sinh(x). 

3. If x ^ 0.481212, then use sinh(x) = % (V - ~V). 

Appendix F: Algorithms 55 



4. If fg x < 0.481212, then compute sinh(x) by use of a polynomial: 

sinh(x) 

= a + «i x 2 + #2* 4 + 03X 6 + aix 8 + a 5 x 10 . 

The coefficients were obtained by the minimax approximation (in relative 
error) of sinh(x)/x over the range ^ x 2 ^ 0.23156 under the constraint 
that a shall be exactly 1.0. The relative error of this approximation is less 
than 2- 5607 . 

5. The value of cosh(x) is computed as cosh (x) = sinh|x| H — ^. 

e 

This computation uses the real exponential subprogram (dexp) and the real 
sine/ cosine subprogram (dsin/dcos). 

Effect of an Argument Error 

To understand the effect of an argument error upon the accuracy of the answer, 
the programmer must understand the effect of an argument error in the dsin/dcos, 
dexp, and dsinh/dcosh subprograms. 



Square Root Subprograms 
SQRT 

Algorithm 

1. If x = 0, then the answer is 0. 

2. Write x = 16 2p_(7 • m, where 2p — q is the exponent and q equals either or 1; 

in is the mantissa and is within the range — ^ m < 1. 

16 

3. Then, y x = 16" • 4"«V™. 

4. For the first approximation of V x, compute the following: 



/ 1.288973 \ 

* = !*' 4- -(1.681585 ____). 



This approximation attains the minimax relative error for hyperbolic fits of Y"x. 
The maximum relative error is 2~ 5 - 748 . 
5. Apply the Newton-Raphson iteration 



•' /n+i = 1/2 ( y - + i) 

twice. The second iteration is performed as 

(x \ x 
Vi 1+ — , 



with a partial rounding. The maximum relative error of y 2 is theoretically 

2-25.0. 

Effect of an Argument Error 
1 

DSQRT 

Algorithm 

1. If x = 0, then the answer is 0. 

2. Write x = 16 2p -* • m, where 2p — q is the exponent and q equals either or 1; 

m is the mantissa and is within the range — g m < 1. 

16 
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3. Then, y/x = 16* • 4-« V^. 

4. For the first approximation of V *> compute the following: 

y Q = 16* • 4 1 -* • 0.2202 (ro + 0.2587). 
The extrema of relative errors of this approximation f or q = are 2~ 3 - 202 at 

m = 1, 2- 8 - 285 at m = 0.2587, and 2~ 2 ° 25 at m = — . This approximation, rather 

lo 

x 
than the minimax approximation, was chosen so that the quantity j/ 3 be- 

low becomes less than 16* -8 in magnitude. This arrangement allows us to 
substitute short form counterparts for some of the long form instructions in the 
final iteration. 

5. Apply the Newton Raphson iteration 



lfr + t -%(*+ y) 



four times to t/ , twice in the short form and twice in the long form. The final 
step is performed as 



Vi = «/3 + Vz 



(£-») 



with an appropriate truncation maneuver to obtain a virtual rounding. The 
maximum relative error of the final result is theoretically 2~ e3 - 23 . 

Effect of an Argument Error 

1 

£ ~T 8 

CSQRT/CDSQRT 

Algorithm 



1. Write V x + iy = a + i&. 

I i^i + i^ -j- {|#i 

2. Compute the value % = *| r as k • V ti>i + w 2 where k y Wt and w 2 

are defined in 3, or 4, below. In any case let v x = max ( |x|, |i/j ) and 

©2 = min (|x|, |^|). 

3. In the special case when either v 2 = or v t > o 2 , let u>i = t? 2 and tv 2 = t> x so 
that Wx + u> 2 is effectively equal to v t . 

Also let k = 1 if v x = |x| and 

ft = 1/V2tft>! = |y|. 

4. In the general case, compute F = ^jVi -^ Vii— I . 

If |x| is near the underflow threshold, then take 

o>i = |x|, u? 2 = t>i • 2F, and fc = 1/ Y2. 
If v x • F is near the overflow threshold, then take 

itfi = |x|/4, u; 2 = ©i -F/2, and fc = y/2. 

In all other cases, take u>i = |x|/2, u? 2 = v x • F, and k = 1. 

5. If * = 0, then a = and b = 0. 

If z ^ and x ^ 0, then a = z, and 

2z 



If z =^ and x < 0, then a = 



, and 
• z. 
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2z 
fc = (sign!/) «z. 



The algorithms for both complex square root subprograms are identical. 
Each subprogram uses the appropriate real square root subprogram (sqrt or 

DSQRT ) . 

Effect of an Argument Error 

The effect of an argument error depends upon the accuracy of the individual 
parts of the argument. If x + iy = r • & h and V x + iy = R • e iH , 

thenc(-R) ~ -y 8(0,and c (ff ) ^ 8(h). 



Tangent and Cotangent Subprograms 

TAN/COTAN 

Algorithm 

TT 

1. Divide |x| by — and separate the result into integer part (q) and the fraction 

TT 

part (r). Then \x\ = — (q + r). 

2. Obtain the reduced argument (to) as follows: 

if q is even, then to = r 

if q is odd, then to = 1 — r. 

The range of the reduced argument is ^ w :g 1. 

3. Let q> = g mod 4. 

Then for g = 0, tan \x\ — tan I — • to 1 and cot |xj = cot ( — • w 1, 

g () = 1, tan \x\ = cot I — • w I and cot |x| = tan I — • tv I, 

q — 2, tan |x| = — cot I — • tv J and cot |xj = — tan I — • to I, 

qr = 3, tan |%; = — tan I — • tv J and cot \x\ = — cot I — • tv \. 

t. The value of tan f -^- • u> J and cot f -^- • to J are computed as the ratio of two 
polynomials: 

/ * \ wP(u) (<* \ Q(u) 

tan ( —r- • W I = ~7T7 — r— , COt 1 —7- * IV I ~ =77 — r 

\ 4 / Q(u) ' \ 4 / wP(u) 

wheiv; a = yztv' 1 a nd 

F(u) =■ - 8.460901 + u 

Q( u ) = - 10.772754 + 5.703366 • u - 0.159321 • u 2 . 

These coefficients were obtained by the minimax rational approximation (in 
relative error) of the indicated form. The maximum relative error of this 
approximation is 2~ 26 . Choice of u rather than to 2 as the variable for P and Q 
is to improve the round-off quality of the coefficients. 

5. If x < 0, then tan(x) = — tan jx|, and cot(x) = — cot |x|. 

6. This program is provided with two kinds of error controls. One is for arguments 
whose magnitude is greater than 2 IK • tt. The other is for arguments which are 
very close to a singularity of the function. In either case, the precision of the 
argument is deemed insufficient for obtaining a reliable result. More specifically, 
the second control screens out the following arguments : 

tf) \ x \ = 16~ 63 for cotan (the result would overflow). 

b ) x is such that one can find a singularity within eight units of the last digit 
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value of the floating-point representation of the sum q + r. Singularities are 

cases when the cotangent ratio is to be taken and w = 0. 
The test threshold of this control can be dynamically modified by assembler 
code programs. 

Effect o# an Argumont Error 

A 2 

E ~ tt~ 7 » and c ~ . /rt x for tan(x). Therefore, near the singularities 

cos 2 (x) sm(2x) v ° 

x = I fc +— I 7T, where fc is an integer, no error control can be maintained. This 
is also true for cotan(x) for x near kw, where k is an integer. 

DTAN/DCOTAN 

Algorithm 

1. Divide |x| by ~ and separate the result into integer part (q) and the. fraction 

part (r). Then |x| = -J(g + r). 

2. Obtain the reduced argument ( to ) as follows ; 

if q is even, then w = r 
if g is odd, then to = 1 — r. 

The range of the reduced argument is ^ w :§ 1. 

3. Let q$ = q mod 4. 

Then for (70 = 0, tan |x| = tan ( — • to land cot \x\ = cot I — • to 1, 
g = 1, tan |x) = cot I — • to land cot |x| = tan 1 — • to 1, 
g = 2, tan |x| = — cot ( — • tv land cot |x| = - tan ( -j- • w I, 
go = 3, tan |x| = - tan (~ • to Jandcot |x| = - cot ( JL . w 1 

4. The value of tan 1 — • to \ and cot I -r- • ti; | are computed as the ratio of 
two polynomials: 

tan ^_.^^___ andcot J_.«,J a ^^y^y. 

where both P and Q are polynomials of degree 3 in to 2 . The coefficients of P 

and Q were obtained by the minimax rational approximation (in relative error) 

1 / IT \ 

of — tan 1 — w 1 of the indicated form. The maximum relative error of this 

approximation is 2~ 55 - 6 . 

5. If x < 0, then tan(x) = - tan |x|, and cot(x) = - cot |x|. 

6. This program is provided with two kinds of error controls. One is for argu- 
ments whose magnitude is greater than 2 50 • it. The other is for arguments which 
are very close to a singularity of the function. In either case, the precision of 
the argument is deemed insufficient for obtaining a reliable result. More 
specifically, the second control screens out the following arguments: 

a ) |x| ^ 16- 63 for cotan ( the result would overflow ) . 

b ) x is such that one can find a singularity within eight units of the last digit 
value of the floating-point representation of the sum q + r . Singularities are 
cases when the cotangent ratio is to be taken and to = 0. 
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The test threshold of this control can be dynamically modified by assembler 
code programs. 

Effect of an Argument Error 

A '- 2 

E ~ ^t-t, and e ~ . , n s for tan(x). Therefore, near the singularities of 

cos 2 (>) sm(2x) v ' & 

x = I it -i — — I ?r, where k is an integer, no error control can be maintained. 

This is also true for cotan(x) for values of x near hrr, where k is an integer. 



Implicitly Called Subprograms 

The entry point names of the following implicitly called subprograms are gener- 
ated by the compiler. 



Complex Multiply and Divide Subprograms 

CDVD#/CMPY# (Divide/MultipJy for COMPLEX*8 Arguments) 
CDDVD#/CDMPY# (Divide/Multiply for COMPLEX*16 Arguments) 

Algorithm 

Multiply: (A + Bi) (C + Di) = (AC - BD) + (AD + BC)i 
Divide: (A + Bi) I (C + Di) 
1. If |C[g |D|, set 

A = B,B = - A,C = D 7 D = - C, since 

A + Bi B = Ai 



C + Di D - a 



before step 2. 



A B D 

2.SetA^-B^-D^- ; 

then compute 

A + Bi A f + B'i A' + B'D B' - A f D' 



vrv. "*" 



C + Di 1 + D'i 1 + D f D f - 1 + D'D f 

Error Conditions 

Partial underflows can occur in preparing the answer. 

Complex Exponentiation Subprograms 

FCDXI# (COMPLEX* 16 Arguments) 
FCXPi# (COMPLEX*8 Arguments) 

Algorithm 

The value of y x -f y4 = (z x + z 2 i) j is computed as follows. 

K 

Let i/l = £ r k * 2 k where r /i: = or 1 for fc = 0, 1, . . ., K. 
k = o 

Then z | j i = 7r z 2 \ and the factors z 2k can be obtained by successive squaring. 

r k * 

More specifically: 

1. Initially; k = 0, n< (,> = |/|, t/^°> + t/ 2 <°>i = 1 + Oi, 

Zi (o } + Z2 (o )f = Zl + Z2 i 
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2. Raise the index fc by 1, and let n**" 1 * = 2q + r, where q is the integer 
quotient and r = or 1. 

3. Let n<*> = <?. 

4. If r = 0,thenffi<*> + */2 (fc) i = id**- 1 * + Ifc**- 1 '*. 

Iff = l,thenjfi<*> + y 2 (fc) i =(t/i (fe ~ 1) + Ife**" 1 **) (ssi** -1 * + afe'*- 1 **)- 

5. If n<*> +0,thenz l <*> + z 2 (fc) * = (z^ k -^ + 2 2 (fe " 1) ) 2 s and steps 2 
through 5 are repeated until n ik) = 0. 

6. Whenn<*> = 0, and /> 0, then y 2 + t^i = yi ik) +y 2 (fc) f- 
If; < O^hent/x + y 8 < = (1+ W) / (yi {k) + t/ 2 w 0. 



Exponential /on of a Reaf Base to a Real Power Subprograms 

FDXPD# (REAL*8 Arguments) 
FRXPR# (REALM Arguments) 

Algorithm 

1. If a = and b ^ 0, error return. 

If a = and fo > 0, the answer is 0. 

2. If a 7^ and fo = 0, the answer is 1. 

3. All other cases, compute a & *as e & • loff a . In this computation the exponential sub- 
routine and the natural logarithm subroutine are used. If a is negative or if 
b • log a is too large, an error return is given by one of these subroutines. 

Error estimate 

The relative error of the answer can be expressed as (ci + c 2 ) b • log (a) + c 3 
where e u e 2 , and e 3 are relative errors of the logarithmic routine, machine multi- 
plication, and the exponential routine, respectively. 

For fdxpd#, €i ^ 3.5*10- 16 , h <i 2.2xl0- 16 , and €a ^ 2.0xl0- 16 . Hence the 
relative error :§ 5.7xl0~ 16 x | fo • log a | + 2.0xl0 -16 . Note that b • log a is the 
natural logarithm of the answer. 

For frxpr* €i ^ 8.3xl0- 7 , € 2 ^ 9.5xl0~ 7 , and e 3 ^ 4.7xl0~ 7 . Hence the relative 
error g 1.8x10" fl x | b • log a | + 4.7xl0~ 7 . 

Effect of an Argument Error 

[fl( 1 + Si )] b( 1 + So) ss a h ( 1 + 8 2 b • log a + fcS t ). Note that if the answer does 
not overflow, \b • log a\ < 175. On the other hand b can be very large without 
causing an overflow of a b if log a is very small. Thus, if a = 1 and if b is very 
large, then the effect of the perturbation 8 t of a shows very heavily in the relative 
error of the answer. 



Exponentiation of a Real Base to an Integer Power Subprograms 

FDXPI# (REAL*8 Arguments) 

Algorithm 

K 

The value of y = a' is computed as follows: Let |/| = ]T r k 2 k where r k = or 1 

k = 

for fc = 0, 1, . . ., K. Then al-'l = tt a 2 * and the factors a- k can be obtained by suc- 

cessive squaring. 
More specifically: 
1. Initially: k = 0, n<°' = |/|, y m = l,-andz<°> = a. 
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2. Raise the index k by 1, and decompose n ik ~ u = 2q -f r, where g is the 
integer quotient and r = or 1. 

3. Let n ik) = q. 

4. If r = 0, then !/<*'» = y u ~ l \ 

If r = l,thent/ a> = ^/.--i^a-n. 

5. If n ik) 4= 0, then z ik) = z ( *~ 1) z u ' 1 " ) , and steps 2 through 5 are repeated 
until n ,A> = 0. 

6. When ?i'*-» = 0, and / ^ 0, then y = y ik K If / < 0, then y = -j^. 

y 
Note : The negative exponent is computed by taking the reoiprocal of the posi- 
tive power. Thus it is not possible to compute 16.0** — 64 because there is a lack 
of symmetry for real floating-point numbers — i.e., 16.0** — 64 can be represented, 
but 16.0**64 cannot. The result is obtained by successive multiplications and is 
exact only if the answer contains less than 14 significant hexadecimal digits. 

FRXPI# (REALM Arguments) 
Algorithm 

Tins subprogram has the same algorithm as fixpi#, which follows. 



Exponentiation of Integer Base to Integer Power Subprogram 

FIXP1# (INTEGERM Arguments) 

Algorithm 

The value of L = V is computed as follows: Let / = Y r k • 2 k where r k = 

k = 

or 1 for k : = 0, 1, . . ., K. Then V = tt P\ and the factors P k can be obtained by 

successive squaring. 
More specifically: 

1. Initially: k = 0, n MM = /, y ,,M = 1, and m ilu = 1. 

2. Raise the index k by 1, and decompose n a * -1> = 2q + r, where q is the 
integer quotient and r = or 1. 

3. Letn a > = q. 

4. If r - 0, theni/ a » = tj a - ,( . 



If r = l,theni/ a » = y ik 



-n . m i*--n 



5. If n Un -.-\= 0, then m iki = m a '- u • m u "" ], 5 and steps 2 through 5 are repeated 
until ti u ' : = 0. 

6. Whenn a » = 0, L = L u '». 

iVoto: The result is obtained by successive multiplications. The result is exact 
only if it is less than (2**31) — 1. Results are meaningless when this limit is 
exceeded and may even be of changed sign. 
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